Neutrinos with Lorentz-violating operators of arbitrary dimension 
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The behavior of fermions in the presence of Lorentz and CPT violation is studied. Allowing for 
operators of any mass dimension, we classify all Lorentz-violating terms in the quadratic Lagrange 
density for free fermions. The result is adapted to obtain the effective hamiltonian describing the 
propagation and mixing of three flavors of left-handed neutrinos in the presence of Lorentz violation 
involving operators of arbitrary mass dimension. A characterization of the neutrino coefhcients for 
Lorentz violation is provided via a decomposition using spin-weighted spherical harmonics. The 
restriction of the general theory to various special cases is discussed, including among others the 
renormalizable limit, the massless scenario, flavor- blind and oscillation- free models, the diagonaliz- 
able case, and several isotropic limits. The formalism is combined with existing data on neutrino 
oscillations and kinematics to extract a variety of measures of coefficients for Lorentz and CPT vi- 
olation. For oscillations, we use results from the short-baseline experiments LSND and MiniBooNE 
to obtain explicit sensitivities to effects from flavor-mixing Lorentz-violating operators up to mass 
dimension 10, and we present methods to analyze data from long-baseline experiments. For prop- 
agation, we use time-of-flight measurements from the supernova SN1987A and from a variety of 
experiments including MINOS and OPERA to constrain oscillation-free Lorentz-violating operators 
up to mass dimension 10, and we discuss constraints from threshold effects in meson decays and 
Cerenkov emission. 



I. INTRODUCTION 

Nongravitational phenomena are well described by the 
minimal Standard Model (SM) of particle physics [1]. 
However, the existence of physics beyond the SM is es- 
tablished by the confirmed observation of neutrino oscil- 
lations, while the SM itself is believed to be a low-energy 
effective theory emerging from an underlying unified de- 
scription of gravity and quantum physics at the Planck 
scale. Potential experimental signals exposing founda- 
tional Planck-scale physics are therefore of great interest 
but are challenging to identify. One proposed class of 
signals involves the breaking of Lorentz symmetry asso- 
ciated with tiny deviations from relativity [2]. In recent 
years, searches for Lorentz violation and related CPT vi- 
olation have been performed with a wide range of systems 
and at impressive sensitivities [3]. 

Effective field theory can be used to describe Lorentz 
violation in realistic models at attainable energies [4] . In 
this approach, CPT violation is associated with Lorentz 
violation [5]. The comprehensive effective field the- 
ory incorporating the SM and General Relativity and 
characterizing general Lorentz violation is the Standard- 
Model Extension (SME) [6, 7]. Its Lagrange density 
behaves as a scalar density under observer transforma- 
tions, so the SME action is coordinate independent. 
Each Lorentz-violating term is formed by contracting a 
Lorentz-violating operator of a given mass dimension d 
with a controlling coefficient. The mass dimension d pro- 
vides a partial classification of the operators, and it offers 
a rough sense of the size of associated experimental ef- 
fects. 

In this work, we focus on Lorentz violation in neutri- 
nos. The interferometric nature of neutrinos and their 



tiny apparent mass scale makes them natural probes 
for Planck-scale effects such as Lorentz and CPT vio- 
lation. Various searches for Lorentz violation with neu- 
trinos have achieved sensitivities to SME coefficients at 
levels comparable to Planck-suppressed effects, including 
ones by the LSND [8], Super-Kamiokande (SK) [9], ML 
NOS [10, 11], IceCube [12], and MiniBooNE [13] collabo- 
rations, while evidence for superluminal neutrinos has re- 
cently been presented by the OPERA collaboration [14] . 

The introduction of neutrino coefficients for Lorentz 
and CPT violation [6] has led to numerous efforts 
to understand their implications for neutrino behav- 
ior. The general effective hamiltonian describing neu- 
trino propagation and mixing in the presence of Lorentz- 
violating operators of renormalizable dimension contains 
four types of coefficients [15], leading to many novel ef- 
fects that can be revealed in suitable experiments [16, 17]. 
One interesting theoretical challenge is the construction 
of a global Lorentz-violating model describing all estab- 
lished neutrino behavior and perhaps also one or more 
of the known anomalies, without the two usual neutrino 
masses [18-22]. Many of the SME-based phenomenologi- 
cal studies of neutrino behavior focus on the special case 
of isotropic Lorentz violation in a preferred frame [23-44] , 
sometimes in the two-flavor limit. Several works treat 
anisotropic effects without lepton-number violation [45- 
52]. Models also exist that incorporate nonconservation 
of lepton number, which leads to neutrino-antineutrino 
mixing [53-55]. The propagation of free neutrinos in the 
SME follows geodesies in a pseudo-Riemann-Finsler ge- 
ometry [56]. 

The primary goal of this paper is to extend available 
techniques for handling Lorentz violation in the neu- 
trino sector to include operators of arbitrary mass di- 
mension d. We are motivated partly by the notion that 
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the usual SM represents the dominant component of a 
low-energy effective theory, with subdominant terms in- 
volving Lorentz violation given by the minimal SME with 
d = 3 OT d = 4. In this picture, higher-order corrections 
involve Lorentz- violating operators of larger d, which are 
expected to grow in significance as energies increase and 
may thereby provide a link to the underlying theory [57]. 
In some cases, such as supersymmetric Lorentz- violating 
theories [58] , operators with d > 4 can represent the dom- 
inant corrections to the SM. Similarly, in noncommuta- 
tive quantum electrodynamics [59] , the action written in 
terms of the conventional fermion and photon fields con- 
tains only SME operators of dimension d > 6 [60] . 

The analysis performed in this work has many parallels 
to the treatment for Lorcntz-violating operators at arbi- 
trary d performed for the photon sector [61], and much 
of the methodology established for that work is appli- 
cable here. The formalism contains as a limiting case 
the earlier systematic investigation of neutrino operators 
of renormalizablc dimension [15] and the corresponding 
SME-based studies of neutrino Lorentz violation men- 
tioned above. However, the present work generalizes 
the c;xisting treatment to include effects at leading or- 
der in both mass and Lorentz violation, uncovering novel 
Lorentz- violating effects involving neutrino helicity flip. 
It also incorporates as limiting cases studies of nonrenor- 
malizable neutrino operators [15, 62-68], including dis- 
cussions of modifled dispersion relations in the context 
of super luminal neutrinos [69-99]. 

To achieve a reasonable scope while covering neu- 
trino propagation and mixing in the presence of general 
Lorentz and CPT violation, we take the Lagrange den- 
sity of interest to be quadratic in free fermion fields. Our 
methods and results are applicable to any type of fermion 
and hence are also relevant for other sectors of the SME. 
Implications of this and of possible interaction terms are 
considered elsewhere [100]. In this work, our primary 
focus is the application to multiple generations of left- 
handed neutrinos. We obtain the effective hamiltonian 
describing neutrino propagation and mixing in the pres- 
ence of Lorentz- and CPT-violating operators of arbi- 
trary mass dimension. Since rotations form a subgroup 
of the Lorentz group, careful use of rotational proper- 
ties often benefits experimental analyses, and so we de- 
velop a decomposition in spherical harmonics to charac- 
terize effects. We find that all Lorentz-violating features 
of neutrino propagation and mixing are determined by 
four sets of effective spherical coefficients for Lorentz vi- 
olation, which can physically be distinguished by their 
Dirac or Majorana nature and by their CPT properties. 
Using this classification scheme, various limiting mod- 
els can readily be identified and studied, and explicit 
measurements of coefficients for Lorentz violation can 
be extracted from observational and experimental data. 
Here, we use existing results from experiments on oscilla- 
tions, times of flight, thresholds, and Cerenkov emission 
to tabulate measurements and maximal attained sensitiv- 
ities for Lorentz-violating operators of mass dimensions 



d < 10. Many of our results represent the first avail- 
able constraints on the corresponding Lorentz-violating 
effects. 

The structure of this paper is as follows. The gen- 
eral quadratic action for a set of fermion fields in the 
presence of arbitrary Lorentz and CPT violation is pre- 
sented in Sec. II. The specialization of this analysis to the 
neutrino sector, which allows for multiple flavors of left- 
handed neutrinos, is derived in Sec. III. The spherical 
decomposition of the effective hamiltonian for neutrino 
propagation is considered in Sec. IV. We discuss a variety 
of special cases in Sec. V, including the renormalizablc 
limit, the massless scenario, fiavor-blind and single-fiavor 
models, the diagonalizable case, and isotropic models. 
The results are applied in Sees. VI and VII to extract 
numerous limits on coefficients for Lorentz violation from 
existing data on neutrino oscillations and propagation. 

II. FERMIONS 

The construction of a realistic low-energy effective the- 
ory for fermions that is coordinate independent and de- 
scribes general Lorentz violation can be achieved by 
adding appropriate terms to the conventional fermion La- 
grange density. Each additional term is formed by con- 
tracting a coefficient for Lorentz violation with a tensor 
operator, and all possible terms are included [4, 6, 7]. 
In this work we focus on noninteracting fermions, which 
corresponds to restricting Lorentz-violating terms in the 
action to fermion bilinears. The formalism presented in 
this section holds for arbitrary fermions, but our primary 
interest in subsequent sections lies in applications to neu- 
trino physics for which chiral components are physically 
relevant. 

Consider the case of N spinor fields ipa, where a ranges 

over n spinor flavors a = 1,2,..., A''. To allow for Ma- 
jorana couplings in the construction, it is convenient to 
combine the A'' spinors (,/'a together with their charge con- 
jugates = CV'a into a 2A^-dimensional multiplet of 
spinors. 




where A ranges over 27V values. The redundancy in ^ 
implies that it obeys the relationship 

*c?=C*, C=(J J), (2) 

where the 2N x 2N matrix C is defined in terms oi N x N 
blocks in fiavor space. 

In terms of the spinor we can write the general La- 
grange density incorporating Lorentz and CPT violation 
in the form 

S = J Cd'^x, 

£ = ^'^Ah^'id^^^AB - Mab + Qab)'^b +h.c. (3) 
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The first part of this expression generates the usual ki- 
netic term, while the second part involves an arbitrary 
mass matrix Mab ■ The third part contains the Lorentz- 
violating operator Qab, which is a general 4x4 matrix 
in spinor space and a 2N x 2N matrix in flavor space 
that involves derivatives id^. 

Since the effects from Lorentz violation are generi- 
cally expected to be small, possibly arising as Planck- 
suppressed effects, it is reasonable to treat Qab as a per- 
turbative contribution when necessary. This approach is 
adopted in the present work. In principle, the unper- 
turbed theory could be taken as unconventional if de- 
sired. For example, some models of superluminal neu- 
trinos adopt a kinetic term involving the replacement 
«7^a^ ^ ^757^5^ [69 75, 77, 88, 90]. Models of this 
kind can be incorporated in the present formalism with 
a nonperturbative choice of Qab- 

The hermiticity of C implies that the general form of 
Mab can be written 

Mab = niAB + imsABls, (4) 

where m and are hermitian 2N x 2N matrices. The 
relationship (2) implies the conditions 

m = Crri^C, = Crn^C, (5) 

where the transpose acts in flavor space. 

The operator Qab can in general depend on space- 
time position, either in a prescribed way or through dy- 
namical flelds. For instance, explicit Lorentz violation 
occurs when Qab contains a fixed background with non- 
trival Lorentz properties, while spontaneous Lorentz vi- 
olation can arise if Qab involves field variables with dy- 
namics generating tensor vacuum values. An analysis 
incorporating spacctime dependence would be of interest 
but woiild be burdened by theoretical and experimental 
complexities beyond Lorentz violation, so it is advanta- 
geous to focus on operators that conserve energy and 
momentum. This can be assured by requiring the in- 
variance of the action S under spacetime translations, 
which is achieved when Qab is spacetime independent. 
For spontaneous Lorentz breaking, requiring invariance 
under spacetime translations implies neglecting soliton 
solutions, along with any massive or Nambu-Goldstone 
(NG) modes [101]. The latter can be interpreted as the 
photon in Einstein-Maxwell theory [7, 102], as the gravi- 
ton [103], or as a variety of other forces [104]. 

We remark in passing that spacetime independence 
may be a natural feature of a model, or it may be a 
useful approximation describing dominant contributions 
or averaged effects from known or hypothesized forces 
in the vicinity of the Earth. Even Lorentz-invariant in- 
teractions typically generate effective Lorentz violation 
in this way. For example, couplings to a tiny and pre- 
viously unknown Lorentz-invariant inverse-square force 
in the vicinity of the Earth would generate effective 
Lorentz- violating behavior described by the SME. This 



idea is the basis for some of the sharpest sensitivities 
obtained on torsion to date [105] and on some Lorentz- 
invariant effects from quantum gravity [106]. Models for 
new Lorentz-invariant gravitational interactions are vi- 
able only if they are compatible with the numerous ex- 
isting constraints for Lorentz violation [3]. Experimental 
bounds on spacetime-independent Qab obtained in this 
work therefore also constrain Lorentz-invariant models 
involving neutrinos. 

A decomposition of Qab permits the Lorentz- violating 
operators in C to be classified and enumerated. The spin 
part of Qab can be characterized by expanding in the 
basis of 16 Dirac matrices 77, 

Qab = ^Qab7i 

I 

= Sab + iVABlb + VXbI^. + -^.^abI^Ih + YI'ab'^hi^^ 

(6) 

where the 2N x 2A'' derivative-dependent matrix opera- 
tors Qab =^^6 hermitian in flavor space. The derivative 

dependence can be revealed by expressing each Q^ab ^ 
a sum of operators of definite mass dimension d, 

00 

Q'aB = E 2iT^"^ -"^-^f..fc.. . . .Pa._3, (7) 

where = idf_i. Since each has mass dimension 

one, the coefficients Q(^^^"i"2. .ad-3 j^^yg mass dimension 
4 — d. Following the discussion above, these coefficients 
can be taken as spacetime constants. 

A useful refinement of the above decomposition in- 
volves first splitting Qab as 

YpJab - Mab + Qab = T^^bP-^ - Mab , (8) 

in analogy to the usual split in the single- fermion limit of 
the minimal SME [6]. The combination F^^p^ contains 
all operators of even mass dimension, while Mab con- 
tains all those of odd mass dimension. Expanding these 
combinations using Dirac matrices gives 

Tab = l^SAB + c^B^ft^ + dABlalt, 

+^AB + ^lABl^ + \9fB<ynX, 

Mab = mAB + im^ABlb + ffiAB + im^ABlb 

+aABlf, + b'XBj5l^ + ^Hi^^ga^,^. (9) 

The dimensionless operators c^"^, d^"^ are CPT even, 
while the dimensionless operators e^^, Jab^ Sab 
CPT odd. The remaining operators have mass dimen- 
sion one, with ffiAB, f^zAB, being CPT even and 
^AB' ^AB being CPT odd. Note that all the operators in 
Eq. (9) have counterparts in the minimal SME except for 
niAB and fh^AB, which contain only terms of nonrenor- 
malizable dimension. 
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Since appears contracted with in Eq. (8), the 

operators c^ab, e^s! /ab' 5a'b ^'^^ ^^^^ automat- 
ically contracted with p^. It is therefore natural and 
convenient to define the contracted operators 



Cab 



'-AS^'!^' 



<B 

/as 



/abP^^. 



9ab 



gf^P.- (10) 



The CPT properties of these contracted operators 
matches those of their counterparts in the minimal SME. 
Using this definition reveals the relationships 



Sab = SAB 
^AB — ^AB ~ "ab ' 
' AB — iJAB ^AB 



friAB-, Vab 

■^AB - 



f. 



AB - m^AB, 



'^AB 



^ABj 



(11) 



between the expansions (6) and (8). 

We can also take advantage of the property (2) to sep- 
arate operators into Dirac and Majorana pieces. For each 
operator of mass dimension d in Qab, the property (2) 
yields the constraint 



Q = {-if-^CCQ^C-^C, 



(12) 



where the transpose acts in both spinor and flavor spaces. 
This implies the conditions 

fni. ^ [-iYc{f'"'fC, (13) 

where now the transpose acts only in flavor space. Us- 
ing these results, we can write the component operators 
Qab ill terms of four N x N block matrices that can be 
designated as being of Dirac or Majorana type, 



S = 

V = 

V'' = 

> = 



'Sd 


(- 


Sm \ 


Vd 




Vm \ 


-pt 


(- 





(Vj^,)t (-l)'^(V^)^ 



' D 



' M 



(14) 



In these expressions, all the Dirac-like matrices are her- 
mitian in flavor space. Depending on the mass dimen- 
sion, each Majorana matrix operator is either symmetric 
or antisymmetric in flavor space, 



' M 



(-1) 



(-i)'(v^^)'', 



Vm - 



= {-ir+\At,) 



(15) 



Using the designations (14), each component operator 
in the expansions (9) can also be split into four N x N 
block matrices of Dirac or Majorana type obeying the 
conditions (15). 

Many physical features of fermions are most conve- 
niently understood in terms of a hamiltonian formula- 
tion rather than an approach based on the Lagrange 
density (3). The presence of arbitrary Lorentz violation 
and the concomitant higher-order time derivatives com- 
plicates the construction of the hamiltonian. However, 
we can flnd an effective 2A'' x 2A'' hamiltonian Hab that 
correctly describes the physics at leading order in Lorentz 
violation. Starting with the modified Dirac equation 



(p ■ jSab - Mab + QAB)'i'B = 0, 



(16) 



we can multiply on the left by 70 and then define Hab 
by the condition 

{E6ab - Hab)^b = 1o{p-i5ab - Mab + Qab)*b = 0, 

(17) 

where E = pq. We can thereby identify 



Hab = 1o{p-i6ab + Mab 
= {Hq)ab + SHab, 



Qab) 



(18) 



where {Ho)ab = Jo{p-7Sab+Mab) is the usual hamilto- 
nian with conventional energy Eq and 5 Hab = ^"/oQab 
is the Lorentz-violating perturbation. Note that the lat- 
ter term typically depends on E. However, the changes 
to the energy £^0 induced by SHab are pcrturbativc by 
construction, so at leading order 6Hab can be evaluated 
at the conventional energy Eq. The leading-order effec- 
tive hamiltonian can therefore be written as 



Hab = {Hf))AB 

-70 («5ab + iVABlb 

+^Ab7/. + -^AS757,* + \T'AB^t^^)\E^E. 



(19) 



III. NEUTRINOS 



Our primary interest in this work lies in the neutrino 
sector of the SME. Since the observed neutrinos are chi- 
ral fermions, describing their properties in the presence of 
Lorentz and CPT violation rcqiiires projecting the gen- 
eral formalism presented above onto left-handed fields. 
In what follows, we retain all leading-order terms from 
Lorentz violation arising from operators of arbitrary mass 
dimension, including terms linear in neutrino mass. This 
incorporates and extends our earlier analysis for opera- 
tors of renormalizable dimension [15], which treated as 
negligible all terms involving the product of a coefficient 
for Lorentz violation with a neutrino mass. 

To proceed with the analysis, it is convenient to intro- 
duce left- and right-handed mass matrices rriL and mR 
satisfying mn — {itil)^ = m + im^, which combine to 
form M according to 



M = ttilPl + mnPR, 



(20) 
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where = (1 - 75)72 and Pr = (1 + 75)72 are the 
usual chiral projection operators. The components of 
the matrix mR = m\ can be identified with Dirac- or 
Majorana-type masses by separating mn into four Nx N 
submatrices according to 



TUrC 



L D 
D'^ R 



(21) 



Here, R and L are the right- and left-handed Majorana- 
mass matrices, while D is the Dirac-mass matrix. The 
complex matrices R, L, D arc restricted only by the re- 
quirement that R and L arc symmetric. 

In the absence of Lorentz violation, the general equa- 
tion describing massive left-handed fermions is 



p ■ -ytpL - Ltpl - DtpR = 0, 



(22) 



and mixing between left- and right-handed neutrinos van- 
ishes ii D = 0. For nonzero D, a seesaw mechanism is 
usually invoked to suppress left-right mixing [107], based 
on the assumption that R is large. Since the right-handed 
neutrinos obey 



(23) 



a large R implies « —R~^D^'ipj^ . The behavior of 
left-handed neutrinos is therefore well approximated by 
the equation 



0, 



(24) 



where the effective left-handed mass matrix m; is given 
by 

mi=L- DR'^D^. (25) 

Note that mi is symmetric, m; = mf . 

Since experiment shows that propagating neutrinos are 
left-handed and that any right-handed components play 
a negligible role, applying a left-handed projection pro- 
duces an excellent approximation to the physical neutrino 
behavior. In x block form, the relevant projection 
of the hamiltonian H is 



Hr 



Pl 0\rr(PL 

PRr'yo pr 



(26) 



As usual, the projectors imply that this expression can be 
reduced to an operator acting on two-dimensional Weyl 
spinors. Wc introduce a^^ = (cr", cr^), where (T° is the 2x2 
identity matrix and arc the usual three Pauli matrices 
with adjoint matrices ct'^ = {a^,—a^). Denoting by ^ 
the two-component Weyl spinor associated with a four- 
component Dirac spinor Pl?/', the 2A'^-dimensional multi- 
plet ^ in Eq. (1) can be replaced with a 2A''-dimensional 
multiplet $w of the form 



(27) 



where (f)'^ ~ icr'^cj)* . Flavor indices are suppressed in these 
expressions. Similarly, the hamiltonian iJ^ can be re- 
placed with its Weyl counterpart ifw 

In the absence of Lorentz violation, Eq. (24) becomes 

p-a(j)- micjf = 0, (28) 

and the hamiltonian takes the form 

\ ml p ■ a J 
The Lorentz-violating piece 5H in Eq. (18) becomes 



(5ilw = 



where 



Sl = Sm + iVm, V£ = + 



(30) 



(31) 



Note that the preservation of chirality ensures that the 
orthogonal combinations and Sr are absent from i?vv- 
The full hamiltonian 



i?w = (-ffw)o + 



(32) 



can be block diagonalized within a suitable approxima- 
tion. Wc proceed here treating neutrinos as relativistic 
particles, but performing a nonrelativistic diagonaliza- 
tion of H^/\f could also be of interest in certain contexts 
beyond our present scope. These could include, for exam- 
ple, experiments with neutrinos of ultra-low energy such 
as measurements of the beta-decay endpoint, or studies 
of the cosmic neutrino background. 

For the relativistic case, we can block diagonalize Hy/ 
to order mf using the transformation 



U 



7nip-(T 
2p^ 



1— J- 



(33) 



2p2 



8p2 



Consider first the Lorentz-invariant piece (ifw)o- This 
becomes 



(F^)o = U{H^)oU^ 




3- 



p ■ cr(l -I '■ — - 



ill 

2p2 



which acts on the transformed Weyl doublet 



(34) 



(35) 



To find the effective hamiltonian governing the propaga- 
tion of neutrinos, we expand in helicity components 
as 

0' = [^(i,p)e^-P + P*(i,p)e-'^-P]ep, 
{4>'f = [A*{t,p)e-^--P + B{t,p)e--P]i^, (36) 
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where is a normalized negative-hclicity spinor satisfy- 
ing p cr^p = — |p|Cp- The ampUtude A{t,p) is associated 
witli negative-lichcity neutrinos, while B{t,p) is associ- 
ated with positive- helicity antineutrinos. In the Lorentz- 
invariant limit, restricting to the positive-energy part of 
the Schrodinger equation gives 



. d 



(/ieff)o ( ^ 



(37) 



where 



\P\ 



miml 



mlmij 



is the usual effective hamiltonian for ncnitrino and an- 
tineutrino amplitudes in the Lorentz-invariant limit. 

To obtain the Lorentz- violating contribution 6h to the 
effective hamiltonian, we first diagonalize 6H^j\[ as 

6Hi^ = U 6Hw C/^ 

= 6Hw + [5U,5Hw]+0{mf), (39) 



where 



SU 



p ■ (7 f —mi 



2p2 



t 







(40) 



We then obtain 5h by projecting onto the positive-energy 

piece. 



5h 



^p^tj^^fwl^o 



(41) 



An explicit result can be obtained using the identities 

CK^p = -4^^ep = \/2e^, 
^p^<^4p = -^p^^/x^p = (42) 
where the polarization vector can be taken as 

e^ = ^{0;e, + ie2), e^{-p) = e;{p). (43) 

Here, ei and 62 are arbitrary unit vectors chosen so 
that {p, 61,62} form a right-handed orthonormal triad. 
Adopting the correspondence p = r, ei = 9, €2 = 4> 
to the usual spherical-coordinate unit vectors implies the 
spatial part of is the helicity unit vector e+ introduced 
in Appendix A 2 of Ref. [61]. 

Some calculation along the above lines reveals that to 
order 0(m;) the Lorentz- violating piece 6h of the effec- 
tive hamiltonian takes the form 

5h=^( % "If + %A , (44) 

\P\ \-9ls + Hls -als-clsJ 

where conjugation and transposition are flavor-space op- 
erations. For convenience and clarity, this expression 



splits each N x N hamiltonian block into CPT-odd and 
CPT-cven parts, where the notation reflects the CPT 
properties of the corresponding operators in the minimal 
SMB. The CPT-odd parts take the form 

oieflF = PndL - e; -I- 2ie^e^g^ , 

5eflF = iV2p^e^g'^^ + V2e^af, (45) 

while the CPT-even terms are 

Ceflf = p^l^ -m+ 2ie^e*ff/"', 
^eff = iV2p^e.H^^ + V2e^c<^. (46) 

In these expressions, each quantity TJ(^'^ is defined as the 

combination 7^"^ = ^{T'^" + iT^"), and it obeys the 

identity e^e*?^*^ w —p^T^j^/\p\. Here and below, a 

tilde denotes the usual dual with e^'^"^/2. The operators 
independent of the mass matrix mi are defined as 

o-l = a^D + ^D^ 9m+ = 2 (ffivf + mi), 

r^ = r^ + r^, h^^^ = \ [h^^ + ikz) , (47) 

and they obey the hermiticity and symmetry conditions 

o-l = «)^ 9m+ = ^9m+ = [9m+) , 

= {c^)\ HZ^=tllZ^ = -{H';^^)r (48) 
The operators linear in m; are given by 

fhi = ^{fnM + i'm5M)rnl + ^mi{mM + im5M)\ 

cf = \-Slm,i - \m,i{^if , 

e; = I (cm + ifM)m\ + \mi (cm + ifM) ^ , 

r = hrM^ + hMrM^^ 

Hr = \H^^^m\ + \rm{H^^^)\ (49) 
and they satisfy 



mi = m 



ei = el T = (?r)'> Hr = (i^r 



(50) 



Generically, all the above operators depend on the 4- 
momentum. 

The net effective hamiltonian h^ti is the 2A^ x 2A'' ma- 
trix given as the sum of Eqs. (38) and (44), 



/leflF = {hes)o + Sh. 



(51) 



Note that neutrinos and antineutrinos have identical 
mass spectra despite the presence of CPT violation [5]. 
Note also that the mass-induced operators (49) are in 
principle all determined by the mass-independent opera- 
tors (47) once the mass matrix mi is known. However, 
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coefficient 


d CP T CPT 


M«l--«(i-3 

uai-«<i-3 

^ai...o,_3 

ai...ad_3 

/^i/ai ...arf_3 
9l 

filial. ..aj_3 
Hi 


odd - + - 
even + + + 
odd + + + 
even — + — 
even + — — 
odd - - + 


/ii/ai...a<j_3 
„(ii^ai...a<j_3 
uai...a<j_3 
M«l---a<i-3 


even + — — 
odd - - + 
odd - + - 
even + + + 



TABLE I: Properties of neutrino coefficients under discrete 
transformations. For CP and T, each coefficient must also be 
complex conjugated and multiplied by an additional factor of 
(—1)", where n is the number of spatial indices. 

inspection of Eqs. (45) and (46) reveals that the two kinds 
of operators enter hes through different projections with 
Pf^ and e^, and they therefore represent independent ob- 
servable effects. 

The terms in Sh can be classified according to their 
operator dimension d and their properties under discrete 
transformations. Neutrinos maximally break C and P 
symmetry because these transformations reverse chiral- 
ity, so we consider here only the chirality-prescrving op- 
erators CP, T, and CPT. They transform a Weyl spinor 

according to 

CP : (j)^ (f^{t, x) = r]cp(T'^(l>*{t, -x), 

T: 4>^ (f^{t,x)=rij,a'^4>*{-t,x), 
CPT : -i> (t>^^'^{t, x) = r}cPT<l>{-t, -x), (52) 

where the phases rj^p, rjj, arc arbitrary but combine 
to give rfuprj, = —rj^^pr)^. For definitencss, we choose 
rjcp = 1, ?7t = VCPT = i- Table I summarizes the behav- 
ior of the coefficients for Lorentz violation under CP, T, 
and CPT. Both CP and T act to complex-conjugate each 
coefficient and multiply it by a factor of (—1)", where n 
is the number of spatial indices on the coefficient. The 
first six coefficients listed in the table enter 5h in the 
on-diagonal blocks, which control v -ic^ v and i7 -H- 77 
mixing. The remaining four appear in the off-diagonal 
blocks, which are associated with u i^V mixing. Notice 
that each class of coefficients has a unique set of proper- 
ties. 



IV. SPHERICAL DECOMPOSITION 

Many experimental tests of Lorentz invariance rely on 
searching for anisotropics associated with violations of 

rotation symmetry. Searches of this type require knowl- 
edge of the transformation properties of the coefficients 
for Lorentz violation under rotations. In principle, any 
given rotation can be performed on the cartesian coef- 
ficients for Lorentz violation discussed in the previous 



subsection. However, in practice this may require signif- 
icant calculation, while the results can be cumbersome 
and can disguise basic aspects of rotation symmetry. 

An alternative approach involves decomposing the co- 
efficients for Lorentz violation in spherical harmonics. 
This emphasizes the importance of rotations, and it 
ensures comparatively simple properties under rotation 
transformations. It is useful both in searches for viola- 
tions of isotropy and in theoretical treatments of certain 
models, such as those exhibiting isotropy in a preferred 
frame. A decomposition of this type has already been 
used to classify and enumerate photon-sector operators 
of arbitrary mass dimension [61]. 

To perform this decomposition for neutrinos, we ex- 
pand in spherical harmonics the p^-dependent combina- 
tions appearing in the Lorentz- violating piece (44) of the 
hamiltonian 6h. The terms appearing in the diagonal 
blocks of 6h are rotational scalars, so they can be ex- 
panded using the standard spherical harmonics Yj^ = 
oYjm- For example, the contribution involving the coef- 
ficients can be written as 

P.i^LU = E E'-'--\pr YjmiP) (cf (53) 

djmn 

where the indices a, b range over neutrino flavors as be- 
fore. Note that E can be approximated as E \p\ to 
second order in the small mass m;. Since we are inter- 
ested in Lorentz- violating effects to first order in m; , the 
second-order terms can be discarded and i? « |p| can be 
assumed for the expansion of the Lorentz- violating oper- 
ators. This approximation is analogous to that used to 
obtain the vacuum coefficients in the photon sector [61]. 

Using this approximation, the relevant spherical de- 
compositions for the six types of coefficients appearing 
in the diagonal blocks of dh can be written as 





= Ew" 

djm 








= E \pr 

djm 




{c[y , 




djm 






{ei)ab 


= Ew" 

djm 


-'yjm{p) 


(Jdhab 
Jjm' 


a9l'')ab 


= Eipi'- 

djm 


-^Y^m{p) 


\9l )jm' 


ti^Dab 


= Eipi" 

djm 


-'Yjmip) 





The coefficients (ai'^^)"'. and (ci'^^Y. have mass di- 

mension 4 — d, while the derived coefficients (rrif^Y'. , 

(el'''')" . (qj'^^)"'. , and {H^'^^\°'. have mass dimension 
5 — d. All these coefficients are hermitian in flavor space. 
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coefficient 



d 



number 



CP 



T 



CPT 



\ t / jm 



odd, > 3 
even, > 4 
odd, > 5 
even, > 4 
even, > 4 

d = 3 
odd, > 5 



d - 2 > j > 

d - 2 > j > 

d - 3 > j > 

d - 3 > j > 

d - 2 > j > 

j = l 

d - 2 > j > 



9(d-l)2 
9(d- 1)2 
9(d-2)" 
9(d-2)2 
9(d - If 
27 

9(d - If 



i-^y{Hnt 



\ -1^ / jTn 

i'm) ') . 

\ t / 77TI 



(6'm+) j>n 



H 



(d) \ab 



(d)\ab 



(d)\ab 



even, > 4 
odd, > 3 
odd, > 3 
even, > 4 



d - 2 > J > 1 

d - 2 > j > 1 

d - 2 > j > 1 

d - 2 > j > 1 



12d(d - 2) 
6d(d - 2) 

12d(d - 2) 
6d(d - 2) 



(-1) 



(-iy+™ 



(-1) 



j+m+ll (Jd) -^ab 
Ad) xab 

yjyb 

' — m) 

,(d)\ai> 

' j( — m) 



(d) \ab 



3 + m/ 



(-iy+"^((< 



VSff ) jm 



{d)\ a6 



odd, > 3 
d = 2 
even, > 4 

even, > 2 
odd, > 3 



- 1 > i > 
i = i 

d - 1 > j > 

- 1 > j > 1 

- 1 > J > 1 



9d2 
27 

j2 



12(d^ - 1) 
6(d2 - 1) 



i-iy^-^'{{9iXi-^)), 

i-iy^-{{HiS):U) 



(d)\ba 
jm 



i-^y{ci^)t 

-^y^-{{9iXi-^) 
(-iy^-(iHm:U) 



-0 



^(d)-,a(> 

,(d)Na6 
'eff ) im 



H. 



(d)\ab 



TABLE II: Spherical coefficients and their properties under discrete transformations. 



which implies they obey a relation of the form 



(55) 



The off-diagonal blocks of Sh induce mixing between 
neutrino and antineutrino states with opposite helicities, 
with the operators in the upper right block having helic- 
ity — 1 . Functions with integral helicity can be expanded 
in spin- weighted spherical harmonics gYjm, where the 
spin weight s is the negative of the helicity. The defini- 
tions and properties of the spin-weighted functions used 
here can be found in Appendix A of Ref. [61], and the ex- 
pansion procedure parallels that adopted for the photon 
sector. The result of the decomposition can be written 
as 





= E 

djm 




-|-l^m(p) (S'if+)j„i' 


i^Pn^v{H'^+)ab 


= E 

djm 


pI' 


^-V,y,„(p) (Ot. 


V^e^(af)ab 


= E 

djm 


pI' 




\/2e^(cf)„6 


= E 

djm 




'-\^Y,r^{p){cf^r;i. (56) 



The CPT-even coefficients (il^^)"^ and (cj''^)"^ are 



^(d)\ab 



antisymmetric in flavor space, while the CPT-odd coef- 



ficients (.9M+)jm ^""^^ ('^i''^)jm ^'''^ symmetric. The co- 
efficients and have mass dimension 



4 — d, while the mass- induced coefficients (ft;'^'')^^ and 
/ ((i)\ab j^g^^g dimension 5 — d. 

Since the coefficients (54) and (56) contribute to 5h 
only through the combinations flofr, Qiir, 5ofr, and ffeffj 
experiments are sensitive only to the latter. This sug- 
gests that for practical applications it is useful to consider 
instead the expansions 

djm 

Ew" 

djm 

Eipi'-Vi^.>"(p)(5^«^t. 

djm 

I jm 



<jab 



<^b 
Ceff 



9eS 



Zjab 
-"eff 



-2y. 



(P)(ciff)jl, 



E 

djm 



(57) 



The four sets of effective spherical coefficients appearing 
in these expansions are related to the ten sets of funda- 
mental ones (54) and (56) by 



(d)\ab 



^(d)^ab 

(.9ofr 

(^fj(d)-^ab 



jm 

jm 



\ 1^ / jrn \ /jm v * ^ jm 

,(d) Y''' _±_ f^(d+i')\<^b 



(d) \ab 



jm 



{d+l)\ab 



jm 



(58) 



The reader is cautioned that the d superscript on the 
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effective coefficients for Lorcntz violation may differ from 
the dimension of the underlying operator. Indeed, the ef- 
fective coefficients are typically the sum of a fundamental 
coefRcient labeled by d and a mass-induced coefficient la- 
beled by d -I- 1. Since the latter set is null in the case of 
massless neutrinos, most of the mass effects at first order 
can be absorbed into the zeroth-order terms. The excep- 
tions are the coefficients having the maximal j value of 
j = for which the effects arise purely from the mass 

interactions. For example, the coefficients (c|,ff)"^^j and 

(^eff )iL with superscripts d = 2 arise entirely as mass- 
induced violations from operators of mass dimension 3. 

Table II compiles some properties of the fundamental 
and effective spherical coefficients for Lorentz violation. 
The first six rows concern the operators (54) on the diag- 
onal blocks of 5h, while the next four concern the opera- 
tors (56) on the off-diagonal blocks. The final four rows 
present information about the effective spherical coeffi- 
cients (57). For each coefficient, the table provides the 
ranges of d and j, the number of independent real com- 
ponents for the case of iV = 3 neutrino flavors, and the 
CP, T, and CPT transformation properties. 



V. SPECIAL MODELS 

Specific searches may have enhanced ability to detect 
certain types of spherical coefficients. Moreover, the 
general analysis involves many independent components, 
which may make some studies challenging to perform. 
Special models representing limiting cases may there- 
fore be useful for certain experiments and for theoretical 
purposes such as modeling signals potentially associated 
with Lorentz and CPT violation. 

This subsection considers some special limiting cases 
of the general treatment above. Five classes of limit- 
ing models are discussed. We begin with renormaliz- 
able models, in which all operators have mass dimension 
four or less. The second category is massless models, in 
which all deviations of neutrino oscillation and propaga- 
tion from the usual lightlike behavior can be attributed 
solely to Lorentz violation. Another class is the flavor- 
blind and oscillation-free models, in which either mixing 
is absent or only single-flavor neutrino-antineutrino mix- 
ing occurs. The fourth is diagonalizable models, where 
a constant mixing matrix can simultaneously diagonal- 
ize the mass matrix and all Lorentz- and CPT-violation 
contributions to the effective hamiltonian. Finally, we 
consider several kinds of isotropic models, for which a 
preferred frame preserving rotation invariance exists. For 
definiteness, we take N = 3 neutrino flavors throughout. 



A. Renormalizable models 

The limit of the SME in which the only nonzero 
Lorentz- violating operators are of mass dimension d < 4 



spherical 


cartesian 


j 


number 


/ (4)\ 

(o<^') 

\yctf ) jm 
V ett J jjn 


(3)m (4)m „(4)mi'P 
"■L ^ ' ill 

yM+ 


0,1,2 
1 

n 1 o 

0, 1,2 
1 
1,2 
1,2 


81 
27 

O 1 

81 

36 

96 
48 



TABLE IIL Spherical coefficients for renormalizable models. 

is renormalizable to at least one loop [108]. Theoretical 
aspects of renormalizable SME-based models for neutri- 
nos have been extensively studied [6, 7, 15-18, 20, 22-54] 
and several experimental collaborations have measured 
numerous renormalizable coefficients for Lorentz viola- 
tion [3, 8-13]. Here, adopting a flxed inertial frame with 
cartesian coordinates {t,x,y,z), we establish the linear 
combination of cartesian coefficients that corresponds to 
the spherical ones. Throughout this subsection, fiavor 
indices are suppressed for simplicity. 

To date, discussions in the literature have been re- 
stricted to mass-independent operators of renormaliz- 
able dimension. Inspection of Eqs. (54) and (56) reveals 
that the cartesian coefficients controlling these operators 

are a^^'^, c^^'"', 5^^^"''', and H^^J^'^ . They enter the 
Lorentz- violating piece Sh of the effective hamiltonian in 

the combinations (45) and (46). This structure permits a 
match to the established notation introduced in Eq. (14) 
of Ref. [15], for which the coefficients are conventionally 
denoted as (ai,)^, {cl)^"' , g'"''', and H^"'. The correspon- 
dence is immediate, with a^^^'^ = {ul)^, c^^^^" = {cl)'^" , 

g^^r^^^g'^-^andFi^V^''^//''^ 

The inclusion of mass-induced effects introduces novel 
types of neutrino helicity fiip and thereby leads to ad- 
ditional cartesian coefficients for operators of renormal- 
izable dimension. Using Eqs. (54) and (56), these are 

found to be af cf^'\ ef'^^ , gf^'"' , and Hf^^\ As 
can be seen from Eqs. (45) and (46), they contribute 
to the Lorcntz-violating piece 8h of the effective hamil- 
tonian through momentum and polarization projections 
that differ from those for the mass-independent coeffi- 
cients. This implies the mass-induced coefficients gen- 
erate observationally distinct effects, so they offer ad- 
ditional possibilities for model building along with new 
arenas for experimental searches. 

The cartesian coefficients for operators of renormaliz- 
able dimension are related to spherical ones through the 
expansions (54), (56), and (58). Table III lists the cor- 
respondence between these two sets of coefficients, along 
with the allowed j values and the number of independent 
real components for the spherical coefficients. 

To express these relationships explicitly, it is conve- 
nient to introduce the combinations of cartesian unit vec- 
tors 

a;± = a; =F (59) 
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With this notation, the connections between the spherical are given by 
coefficients {ci'^)j^ and the cartesian ones is 



00 



47r(, 



^ xl (a 



(3)3 



2.9, 

„(4)i 



(4)tUN 



4v^;-(4)tji 

3 yi ' 



10 



'4£/'„(3)2 



-(4)t2t\ 



11 



f 4(a(f)^■-ep 



(3)^ 

cff ) 
(3)N 



(4).i 



'Sir Afe ;;(4)tjfe 
15 y; ' 



2-1 



20 



l;;(4)fjj^ 



5 va; sa; ^ 



a 



(3)N 
eff / 22 



'Stt -fc ~(4)tife 
15 •^ + '^+ i/i 



(60) 



The purely mass-induced coefficients (c^g ) . are given 
by 









-1 


(4) 






10 








11 



167r TT 



I ' 

(3)t2 



3 



(61) 



while the coefficients (c^g) are 



^eff^O 



477 (. 



,(4)t* , 1^(4)«^ 



(4) 



,(4)tj 



11 



(4) 

I'-cff ^2-2 

(4) 



(4)0- 



+ C 



27r -J -fe „(4)ife 
15 ' 



,(4)«j 



(4) 



20 



22 



Xi)zz l(4,)jj\ 



15 V'L 

22: Afe „(4)j's 
15 •^ + -^+ 



„(4)J^ 



(62) 



(4^) 



— /47r „(3)j 



10 



877 J3)z 
3 "i ' 



(4^)ii = 

and the coefficients (5^5 ) are 



47r «(3)j 
3 •^+ "i ' 



(63) 



10 



11 



= I 



(4^) 

(4^) 
(4^) 



2-2 



327r (J4:)tzt _ i {i)txy j „(4)tya;\ 
I.I/M+ 2yM+ "T 2^]^+ ^' 

■ /l67 /'„(4)tjf _ 1„(4)<2J 1 (4)<J2N 

'V 3 '''+ \yM+ 2yM+ "T 2yM+ j' 
5 .yM^ 



(5rff02-i 

(4) 



20 



21 



(4)22 ~ 



/247r ('„(4)«zz _ l_(4)tjj\ 
'V 5 ^.yM+ 3i'M+ ^' 

-i.f^^i i-k „(4)yfc 



(64) 



Lastly, the coefficients {H^g)-^ are given by 
(^i^^)i-i = 



10 



jm 



2W 



(3)^ 

eff /II 



47r -.J 



(2ziJ, 



(3)tj 
M+ 



2-2 



■n ^3 Afe „(4)j*: 
5 ' 



(4)it lAi)z3 , 1 J4)jzN 



(^eff)2-l = ~ 

(iJ^^h - - 
>--"effJ20 



5 



^(4)zj 



2^1 



.(4)i2 



2''i 



6^/^4)22 l„(4)jj^ 



H. 



(3)N 
eff /22 



y^4(cW-+cp^), 



ir ^A: „(4)ife 



Massless models 



(65) 



The ofF-diagonal spherical coefficients (Sog)^^, 
(sitt) jm- ^^'^ (^cff )jm control neutrino- antincutrino 



(2)\ 



mixing. The purely mass-induced coefficients [g^g 



■ jm 



An interesting case of the general formalism is the 
massless limit, to; — >■ 0. Given the compelling experimen- 
tal evidence that neutrinos oscillate amassed in recent 
years, it is reasonable and conservative to adopt the per- 
spective that these oscillations arise from a small nonzero 
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coefficient 


d 


j 


number 


\ ^ / jm 

(5V+)jm 
I Tr{d) \ab 


odd, > 3 
even, > 4 
even, > 4 
odd, > 3 


d - 2 > j > 
d - 2 > j > 
d - 2 > j > 1 
d - 2 > j > 1 


9(d- 1)2 
9(d- 1)2 
12d(d - 2) 
6d(d - 2) 



TABLE IV: Spherical coefficients for massless models. 



neutrino mass matrix. However, the frequently encoun- 
tered claim that oscillations prove neutrinos have mass 

is false, as oscillations can arise from Lorentz and CPT 
violation even when all masses vanish. Massless models 
may therefore be of interest for model building. They are 
also relevant in the ultrarelativistic limit, where masses 
can be neglected but Lorentz-violation operators of di- 
mension four or greater remain important. 

When TO; — > 0, the mass- induced operators m;, e;, 
a^, c^, , Hj^" and the associated coefhcients in Eqs. 
(54) and (56) all vanish. The effective spherical coeffi- 
cients appearing in Eq. (58) therefore reduce in mass- 

less models to {a^)^ = {c^)^, {c^)% = 

{9^)t = i^Ml)t {Hi^)t = This 
implies the angular-momentum index labeling the effec- 
tive coefficients is limited to j < d — 2 rather than to 
j < d — 1. Table IV lists the coefficients for massless 
models, together with the allowed range of d and j and 
the number of independent real components. 

Lorentz- violating massless models studied in the liter- 
ature include the bicycle model [15, 18], its generalization 
by Barger, Marfatia, and Whisnant [20], and the isotropic 
subset of the minimal SME [22] . All these massless mod- 
els involve operators of renormalizable dimension, and 
they can reproduce many observed features of neutrino 
oscillations. However, to date no fully satisfactory mass- 
less model has been presented. The primary issue faced 
by model builders is simultaneously reproducing both the 
KamLAND data [109] and the observed shape of the 
solar neutrino spectrum [110] in the energy range 1-20 
MeV. The KamLAND results can be reproduced using 
the massless Lorentz- violating seesaw mechanism [15], 
but the ~1 MeV scale at which this must be triggered 
is challenging to reconcile with the ~10 MeV scale at 
which the solar-neutrino survival probability passes from 
higher to lower values. The isotropic minimal SME can- 
not accommodate both these features [22], whereas even 
a single mass parameter suffices [21]. 

Given the above issues, it is worth emphasizing that 
the existing literature concerns only a tiny portion of 
the available coefficient space for massless models. The 
potential role of direction-dependent coefficients, includ- 
ing the coefficients (511+)^^ and (-ff|^'_,_)^^ for operators 
of renormalizable dimension, remains largely unexplored. 
Also, many new options exist for realistic model building 
using the nonminimal mass-independent operators clas- 



sified in the present work. Constructing a phenomeno- 
logically viable massless model for neutrino oscillations 
remains an interesting and worthwhile open challenge. 

C. Flavor-blind and oscillation- free models 

A particularly simple limit of the general formalism 
is the flavor-blind limit, obtained by assuming that the 
mass-squared matrix and the Lorentz violation affect all 
flavors in the same way. This limit is unrealistic as 
a global description of neutrinos because no neutrino- 
neutrino oscillations appear. However, under suitable 
circumstances it may represent a useful approximation 
to the physics of neutrino propagation. For example, it 
can be physically relevant when the dominant effects of 
Lorentz violation are flavor blind, with oscillations be- 
ing comparatively small. This can arise via numerical 
values of coefficients or under suitable physical circum- 
stances such as ultra-high neutrino energies. A flavor- 
blind treatment may therefore be appropriate for time- 
of-flight studies, for instance. The flavor-blind cases are 
also useful as toy models of Lorentz- violating effects and 
as a stepping stone to the more general models considered 
below. Note that experimental sensitivity to oscillation- 
free effects is generically reduced because no interferome- 
try is involved. In this subsection, we consider two classes 
of flavor-blind models distinguished according to whether 
they allow single-flavor neutrino-antineutrino oscillations 
or are oscillation free. 

1. Flavor-blind and single-flavor models 

The effective hamiltonian for the flavor-blind limit 
is the restriction of Eq. (51) to three copies of a single 
flavor. Since splits into three identical pieces, we can 
suppress the labeling of the three flavors or consider only 
a single flavor. The effective hamiltonian then takes the 
form 

- + 2\p\ \p\ + Ipl 1^-^:,.^. -S. J • (''^ 

The coefficients Hm+, e;. Hi are absent in this limit be- 
cause Cm, /m, Hm+ are antisymmetric in flavor space. 
The effective components ags and CeS are real, while ^^eflf 
is complex. 

Although neutrino-neutrino oscillations are absent 
in flavor-blind models, the coefficient geff generates 
neutrino-antineutrino mixing for each flavor. This is a 
CPT-odd effect. Even if the coefficient g^vrf vanishes, 
the mass-induced coefficient can contribute to Qes and 
induce oscillations between neutrinos and antineutrinos. 
Mixing is absent only in the CPT-even limit and in the 
special CPT-violating case with both g'^+ and m; van- 
ishing. Examples of single-flavor models with only renor- 
malizable coefficients are presented in Sees. IV B 1 and 
IV B 2 of Ref. [15]. 
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coemcient 


a 


j 


number 


y"-ib )jra 
V fb 1 jm 


odd > 3 


d- 


1 > 


j 


> 


d^ 


even, > 4 


d- 


2 > 


j 


> 


(d-lf 




even, > 2 


d- 


1 > 


3 


> 1 


2{d^ - 1) 


('^of )jm 


odd, > 3 


d- 


2 > 


3 


> 


(d-lf 


\ of J jm 


even, > 4 


d- 


2 > 


3 


> 


(d-lf 



TABLE V: Spherical coefficients for flavor-blind models. 



The flavor-blind hamiltonian (66) can be diagonalized 

in the form 



eff 




The eigenvalues are 



\mi 



2|p 



with 



and the components of the mixing matrix are 

eff o 5cff 



c 



\ + a. 
2A 



fleffj 



(67) 



(68) 



(69) 



(70) 



For experimental investigations, it is \iscful to have a 
description of flavor-blind models accounting for proper- 
ties under rotation transformations. As before, this can 
be achieved by decomposition into spherical harmonics. 
The effective components Uge, Ces, %s can be expanded 
in flavor-blind coefficients as 



fleff = 



Ceff = 



ges 



djm 
djm 
djm 



jm' 



(71) 



where 



(-l)'"(Cfb^)j(_r„)- 



j(-m)' 



(72) 



Table V lists the allowed ranges of d and j for these 
coefficients, along with the number of independent real 
components they contain. 



2. Oscillation-free models 

For certain physical applications and to gain intuition 
within a simple theoretical framework, it can be useful 



to restrict attention to coefficients that cause no mix- 
ing at all. These oscillation-free models are achieved 
starting from the above flavor-blind models and impos- 
ing vanishing neutrino-antineutrino mixing. The form of 
the dispersion relation (68) then implies that the gen- 
eral oscillation-free model can be obtained by setting 
to zero the coefficients {gc's)jm.' Oscillation- free mod- 
els therefore amount to flavor-blind models that conserve 
lepton number. The oscillation- free spherical coefficients 

('^o?)jm ^^'^ (''o?)jm ^PP^^'' expansions of the form 
(71) but have index ranges limited to d > 3 and 4 with 
d — 2 > j > 0, as shown in Table V. Note that most of 
these coefficients describe anisotropic effects, so a generic 
oscillation-free model predicts direction dependence and 
sidereal variations in neutrino and antineutrino proper- 
ties. 

Denoting the neutrino energy in oscillation- free models 
by we can expand in spherical coefficients to obtain 



= \p\ + 



\mi\ 



2\p 



Y 



(el?) 



(73) 



djm 



The antineutrino energy E^^ is obtained by changing the 
sign of the coefficients fai?) . . 

One application of oscillation-free models is the study 
of neutrino propagation. A useful concept in this con- 
text is the group velocity = dE°^ /d\p\, which for a 
neutrino becomes 



2p2 



J2{d-3)\pr'Yjr„ip) 

djm 



(74) 



The antineutrino group velocity = dE^^dlp] takes 
the same form but with a sign change for the coefficient 



D. Diagonalizable models 

Another useful simple limit is the class of diagonal- 
izable models, for which the effective hamiltonian h'^g 
is obtained from the general expression (51) by requir- 
ing all terms to be simultaneously diagonalizable. The 
mass-squared matrix is diagonalized using a momenttun- 
independent mixing matrix U, so in the diagonalizable 
limit each Lorentz-violating term must take a special 
form that is also diagonalizable with the constant ma- 
trix U. In the flavor basis, this implies every Lorentz- 
violating operator commutes with all others and with the 
mass-squared matrix. 

The definition of the diagonalizable limit implies that 
in the mass basis the neutrino behavior is governed by 



three copies of the single-flavor limit discussed in the pre- 
vious subsection. The copies are distinct in that they can 
involve different masses and coefficients for Lorentz vio- 
lation. For the effective hamiltonian h'^g in the flavor 
basis, we can therefore write 

.d ^ (U^ \ / 5^ (e^ \ 

\o u^j [s''* cy \0 E^J 

^ \S'^* J \0 U*J ' ^ ' 

where E^, C^, S'^ are all 3 x 3 diagonal matrices obtained 
by combining the three distinct copies of the single-flavor 
results (68) and (70). 

Diagonalizable models offer potentially interesting 
opportunities to construct realistic models for neu- 
trino oscillations and propagation involving perturbative 
Lorentz-violating effects because the conventional mass 
matrix in the three-neutrino massive model (3i/SM) can 
be adopted together with small Lorentz-violating terms. 
It may also be possible to construct more ambitious diag- 
onalizable models in which Lorentz violation plays a key 
role, perhaps completely replacing one or more neutrino 
mass terms in a vein similar to the puma model [21]. Cer- 
tain diagonalizable models may be usciful as toy models 
or as approximations suitable for describing a more com- 
plete theory in specific physical regimes. For example, 
one simple variation arises if U is taken to be the iden- 
tity. Each neutrino in the resulting diagonalizable model 
is then controlled by a different set of single-flavor co- 
efficients and there are no neutrino-neutrino oscillations, 
although CPT-odd neutrino-antineutrino mixing can still 
arise. 

Diagonalizable models with spatial isotropy represent 
a special restriction considered in Sec. VE2 below. The 
more general class of diagonalizable models having non- 
trivial rotation behavior are conspicuously absent from 
the literature. Many interesting signals are predicted by 
these models, such as direction dependence of time-of- 
flight measurements and of oscillations. 

To study the rotation properties of diagonalizable 
models, it is again useful to perform a decomposition in 
spherical harmonics. The treatment proceeds most easily 
by working in the diagonal basis. Denoting indices in this 
basis with primes, the expansion in spherical harmonics 
becomes 













djm 






Ceff 






'-'yMp){cf)t' 




djm 
















djm 







Table VI lists the spherical coefficients for diagonalizable 
models in the diagonal basis. It also provides the range 
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coefRcient 


d 


3 


number 


\ a / jm 


odd, > 3 


d~ 


1 >i 


> 






even, > 4 


d- 


2>i 


> 


■i(d-lf 




even, > 2 


d- 


i>i 


> 1 


6(d" - 1) 



TABLE VI: Spherical coefficients for diagonalizable models. 



of d and j and shows the number of independent real 
components for each coefiicient. 



E. Isotropic models 

The class of isotropic models, sometimes called 'fried- 
chicken' models due to their popularity and simplicity, is 
generated by restricting attention to the comparatively 
few Lorentz-violating operators that maintain rotation 
symmetry. Since observer boosts mix with rotations, any 
isotropic model is well defined only if its preferred ob- 
server inertial frame is specified. All observers boosted 
with respect to this frame see anisotropic effects. A pop- 
ular choice for the preferrcid frame is the frame of the 
cosmic microwave background (CMB), but other choices 
are possible. Note that choosing the CMB frame for any 
isotropic model implies anisotropics in the canonical Sun- 
centered inertial frame [111, 112] and in Earth- based ex- 
periments. 



1. Generic isotropic models 

The expansion (57) of 5h in spherical harmonics is ide- 
ally suited for investigations of generic isotropic models. 
Only coefficients with j = can contribute in the pre- 
ferred frame. Inspection of Table II and Eq. (58) reveals 
that 5ofi' and H^s vaxist vanish, so isotropic models con- 
tain no operators mixing neutrinos with antineutrinos. 
The effective hamiltonian therefore breaks into two 3x3 
blocks, one for neutrinos and one for antineutrinos. 

For neutrinos, we can write 

/i-' = P+ ^TTi + TT ~ TT' (^^) 

2|p| IpI IpI 

where the ring diacritic is used here and below to denote 
isotropic quantities in the preferred frame [15] . The effec- 
tive hamiltonian for antineutrinos is obtained by trans- 
posing in flavor space and changing the sign of the CPT- 
odd terms, 

A, = |p| + !^-^-^. (78) 
The expansion of the Lorentz-violating terms in spherical 
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isotropic type 


coefficient 
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number per d 


generic 


o{d) 

flab 


odd, > 3 


9 




ab 


even, > 4 
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diagonalizable 


o(d) 

a/ 

a' 


odd, > 3 


3 




a' 


even, > 4 


3 


oscillation-free 




odd, > 3 


1 






even, > 4 


1 



plicit choices 



TABLE VII: Spherical coefficients for isotropic models. 



harmonics takes the simple form 



--fib 
Ceff 



d-2 

^ab ■ 



(79) 



The match to the spherical coefficients appearing in the 
expansion (57) is 



"-ab — 


1 




V47r 


'^ab 


1 




V47r 



'00' 



'00' 



(80) 



Both sets of coefficients are hermitian in flavor space, giv- 
ing 9 real degrees of freedom for each value of d. Since the 
CPT-even effects occur only for even d while CPT-odd 

ones occur only for odd rf, only one of the coefficients a'J^ 

and c'^Ji^ is present at any fixed d. Table VII summarizes 
these basic features. 

Despite their simplicity, isotropic models retain suffi- 
cient complexity to offer interesting prospects as global 
models for neutrino behavior. An interesting example is 
the class of puma models [21], which provide viable al- 
ternatives to the 3i^SM as a global description of existing 
neutrino-oscillation data. For instance, the cgdsm puma 
model is isotropic in the Sun-centered frame and is spec- 
ified by three parameters, consisting of one mass and two 
coefficients for Lorentz violation. The ratio of the two co- 
efficients acts like an effective mass at high energies via 
the Lorentz- violating seesaw mechanism [15]. This ra- 
tio and the mass parameter can be chosen to reproduce 
all accepted neutrino oscillation results, while the third 
degree of freedom naturally generates the MiniBooNE 
anomalies [113, 114] that cannot be accommodated in 
the 3z/SM. The effective hamiltonian in the cgasm puma 
model takes the form of Eqs. (77) and (79) with the ex- 



{mim])ab = 



.(5) 

^o5 



'^ab — 




(81) 



and all other coefficients zero. The numerical values giv- 
ing excellent agreement with experimental data are 



= 2.6 X 10-2^ GeV^ 
a(5) = -2.5 X 10^1^ GeV" 
c(^' = 1.0 X 10"^^ GeV"^, 



(82) 



Extensions of this model can also accommodate the 
LSND anomaly [115] and CPT asymmetries of the MI- 
NOS type [116]. Note that the puma models lie outside 
the class of diagonalizable models described in the pre- 
vious subsection because the nontrivial texture (81) in 
flavor space implies the different terms fail to commute. 



2. Isotropic diagonalizable models 

Combining the diagonalizable and isotropic restric- 
tions described in Sees. YD and YE yields a very simple 
class of models. These models must both have simultane- 
ously diagonalizable Lorentz- violating operators and also 
exhibit rotational invariance in a preferred frame. 

For these models, it is convenient to work in the di- 
agonal basis and in the preferred frame, although care is 
required in applications to Earth-based experiments or 
observations because these are boosted relative to any 
preferred inertial frame and therefore necessarily exhibit 
anisotropic effects. In the diagonal basis and preferred 
frame, the energy of a neutrino of species a' can be writ- 
ten in the form 



Eu,a' = \p\ + 



2|P| 



■Elf 



iail^ - c^^)- (83) 



Three coefficients for Lorentz violation appear for each 
d, one for each neutrino species. 

In isotropic diagonalizable models, the same three co- 
efficients control the behavior of antineutrinos at each d, 
but the antineutrino dispersion relation is CPT conju- 
gated. This changes the sign of the a^? coefficients, so 
that the antineutrino energy is 



Ev,a' — 



\P\ 



2|P| 



(84) 
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Since the Lorentz violation in these simple models 
leads to a power series in positive powers of momentum 
|p| with a single coefficient at each d for each species, the 
group velocity Va' = dE,y,a' /d\p\ for each species a' can 
be obtained immediately. For neutrinos, we find 




For antincutrinos, the group velocity = dE-,a' /d\p\ 
takes the CPT-conjugate form with an opposite sign for 
the coefficient a^f\ We remark that provided the co- 
efficients a|jf'' dominate the deviations from lightspeed, 
this CPT-conjugation property implies that neutrinos 
can be superluminal while antincutrinos are sublumi- 
nal or vice versa. This may be useful in attempts to 
model time-of-flight measurements from laboratory neu- 
trinos and supernova antineutrinos. Notice also that al- 
though isotropic violations with operators of mass dimen- 
sion d = 3 alter the phase velocity, these have no effect on 
the group velocity because they generate only a constant 
shift in the neutrino energy. 

3. Isotropic oscillation-free models 

The above construction for isotropic diagonalizable 
models also incorporates another special class of simple 
models, consisting of the isotropic limit of the oscillation- 
free models described in Sec. VC2. In this case, all 
neutrino species are assumed to have the same isotropic 
properties in the preferred frame. 

For this heavily restricted limit, the a' index appear- 
ing in isotropic diagonalizable models can be disregarded, 
and only one coefficient for Lorentz violation appears for 
each value of d. For example, the neutrino energy in 
these isotropic oscillation-free models becomes 

2|P| Y 
and the corresponding group velocity is 

* = 1 - + - 3)l?>r-' i^^'^ - c^'^ (87) 

Note that these isotropic oscillation-frcc models coincide 
with the isotropic flavor-blind models and with isotropic 
single-flavor models because the isotropic requirement 
forces all Lorentz- violating Major ana couplings to van- 
ish. 

Despite the simultaneous conditions of isotropy, di- 

agonalizability, flavor independence, and no oscillations, 
these models can still exhibit CPT violation because all 
their odd-d Lorentz- violating operators are CPT odd. 
Indeed, the even powers of the momentum appearing 
in the dispersion relation for any isotropic flavor-blind 



model are associated with CPT violation and so the cor- 
responding terms for antineutrinos must change sign, a 
potentially important feature for phenomenology that is 
often overlooked in the literature. 

We remark in passing that any of the above simple 
isotropic deformations of the usual dispersion relations 
for neutrinos can lead to physical and observable efi^ects 
only if some other sector is conventional or exhibits differ- 
ent Lorentz violation. Attempts to invoke a common de- 
formation of Lorentz symmetry across all species merely 
generate conventional physics in an unconventional guise 
[6, 7, 61, 117]. Deformed Lorentz transformations that 
depend on different species are discussed in Ref. [61] and 
are naturally described within the SME framework. 

VI. APPLICATIONS TO OSCILLATIONS 

Neutrino oscillations offer a powerful tool for investi- 
gations of physics beyond the SM because their inter- 
ferometric nature makes them highly sensitive to uncon- 
ventional couplings. Next, we apply the formalism devel- 
oped in the previous sections to explore oscillation effects 
due to Lorentz-violating operators of arbitrary dimen- 
sion, and we obtain explicit constraints on a variety of 
coefficients for Lorentz violation for d < 10. 

Any diagonal terms in the effective hamiltonian (51) 
have no effect on oscillations. We can therefore drop 
the diagonal momentum term in {hes)o. Inspecting Eqs. 
(38) and (44) shows that the effective hamiltonian hose 
controlling oscillations is given by 

^ ^ J_ I \'rn\rni + aeff - Cefi H^s - %n \ 
\P\ \ - gig \m\mi - a^g -c^s)' 

(88) 

The oscillation amplitudes can be found from the time- 
evolution operator S{t) = exp(— i/iosc^)- For practical 
applications, the time t can be identified with the exper- 
imental baseline L, so we can write S{L) = cxp{—ihoscL) . 

An exact treatment of oscillations is typically infeasi- 
ble. One potential issue is that possible decay processes 
including neutrino splitting can introduce nonlinear ef- 
fects in certain regimes. However, no nonlinear effects 
from neutrino decay or other processes have been de- 
tected in experiments to date. A linear treatment us- 
ing ftosc is therefore a realistic and feasible approach for 
obtaining robust and conservative constraints. In what 
follows, two regimes of practical interest are considered. 
The first is the short-baseline approximation, which ap- 
plies when the baseline L is short compared to the effec- 
tive hamiltonian ft-osc and so the transition amplitudes are 
small. This approximation is discussed in Sec. VIA. For 
short baselines, S{L) can be expanded in powers of /lose- 
The second regime, considered in Sec. VI B, is the limit 
of perturbative Lorentz violation. It applies when oscilla- 
tions are primarily due to the mass matrix mjm;. In this 
case, Lorentz and CPT violation can be treated as a per- 
turbation on mass-induced oscillations. The illustrative 
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example of two-flavor maximal mixing with Lorentz vi- 
olation, which also offers intuition about CPT breaking, 
is presented in Sec. VIC. 



where a, /3, and 7 arc Euler angles relating the two 
frames. Denoting the sidereal rotation frequency as 
and the local sidereal time as Tq, we obtain 



A. Short-beiseline approximation 

Expanding S{L) in powers of /Iqsc-^ yields simple 
leading-order approximations for the transition probabil- 
ities from flavor a to flavor b a, 



{\m\mi + aeflF - Ceff)^j 



ab 



2L2 

9 ' 



P 



,2 ' 



p_ 



P _ 



(HeS + 9es) 



ab 



2L2 

9 ' 

2 1.2 



P 



,2 ■ 



(89) 



The survival probabilities can be found by summing over 
possible transitions. These equations generalize Eq. (2) 
of Rcf. [16], where attention was restricted to operators of 
renormalizable dimension and to situations where mass- 
induced oscillations are negligible. 

The above expressions can be used to search for the 
unconventional energy and direction dependences asso- 
ciated with Lorentz and CPT violation. By convention, 
the standard inertial frame used to express and compare 
results for the coefficients for Lorentz violation is a Sun- 
centered frame in which the Z axis is aligned with the 
Earth's rotation axis and the X axis points towards the 
vernal equinox [3, 111, 112]. A beam of neutrinos gener- 
ated on the Earth rotates about the Z axis of the Sun- 
centered frame once each sidereal day, so direct analysis 
of neutrino oscillations in this frame requires an expres- 
sion for the beam direction p as a function of time. 

A more convenient approach adopts instead a stan- 
dard laboratory frame in which the x axis points south, 
y points east, and z points vertically upwards [112]. Typ- 
ically, the source or the detector is chosen as the frame 
origin. In this laboratory frame, the beam direction 
p = pjjjj, is a constant vector, while the coefficients for 
Lorentz violation vary in time instead. 

The spherical-harmonic decomposition developed in 
the previous sections is well suited to analyze this sit- 
uation. Taking Sefi as an example and momentarily sup- 
pressing the flavor indices, we can write 



fleflf 



EH 

djm 



((i)\lab 



%m{p,^{a^); 



(90) 



Neglecting effects from the Earth's boost, the coefB- 
cients (a^g)^^ in the laboratory frame are related by 

a time-dependent rotation to the constant coefficients 



(a['2) . in the Sun-centered frame. This rotation can 

V err / 7yn 



be expressed in terms of Wigner matrices Dl^^,{a, /3, 7), 



jm 



= E«''"'"*^*^w(-x)(«S^w. (91) 



where x is the angle between the Sun-frame Z axis and 
the laboratory- frame z axis, corresponding in the north- 
ern hemisphere to the colatitude of the laboratory. The 

(7) 

quantities c?^^/ are the 'little' Wigner matrices. Explicit 
expressions for the Wigner matrices in the conventions 
used here are given in Eqs. (134)-(136) of Ref. [61]. 

Restoring the flavor indices, we obtain the time- 
dependent expression 



= E 

djmm' 



0) 



(92) 



Each term in this equation depends on the dimension of 
the Lorentz- violating operator through the power of 
on the local sidereal time through a harmonic of the side- 
real frequency, on the direction P\g^ through the spherical 
harmonics and the Wigner matrices, and on the coef- 
ficients for Lorentz violation in the Sun-centered frame. 
The equation thereby determines the linear combinations 
of coefficients of Lorentz violation that can be accessed 
by a given experiment for each flavor transition. 

The dependence on the direction p^^-^ of the beam with 
respect to the Earth can be compactly encoded by defin- 
ing the quantities 



E^ 



^jm'(Plab)C^m(-X) 



which obey the reality conditions 

s-!^jm = (-l)'"'^*-s-^j(-m)- 



(93) 



(94) 



The factors s^fjm = sJ^jmiPxah^X) = sJ^jm.{0,<l),x) vary 
with the laboratory polar angles {6, 4>) and the colatitude 
angle x- Applying the above line of reasoning also to the 
other effective coefficients appearing in /Iqsc then yields 



"eff 



■^b 



-yab 
9eS 



Tjab 
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eflF 



^|p|<i-2^.™.eTe(4^)J^, 

dm 

EW^^""*"M4f^)„> 
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(95) 



dm 
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where 



V Oeff/m 
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(96) 



For completeness, we also define an amplitude for mass 

{A'^]), = Umjmir- (97) 

The transition probabilities can then be written com- 
pactly as 



P 



^^|p|.-3^i™.eTe(_4W)^ 

dm 
dm 



p_ 



dm 



dm 



, (98) 



where 



m 



V / m V "off / m V Ccii / rn 

V'^nii/m V-^acff/m V'^Ceff^OT' 

off ) m 

«1)™ + (42J™- (99) 



As an explicit example, consider the study of Lorentz 
violation by the LSND collaboration [8]. This yielded 
constraints on amplitudes for 77^ — > oscillations with 
sidereal harmonics m < 2. We discuss here the reported 
quadratic bound on the squares of the amplitudes, which 
in current terminology and for fixed dimension d can be 
written as 



|p|2(.-3)A(_4(a)) 



V LSND 

'o 



= 10.5 ± 2.4 ± 1.4 X 10"^^ GeV. (100) 



The LSND analysis assumed that the relevant terms 
in the effective hamiltonian are real, so we must im- 
pose the conditions (aj^fr^)^^^_^^ = (-l)'"(aiff )^|^ and 

(4ff )ir-m) = (-l)'" (ceff ) ■!! • Combining the la errors 



in quadrature gives a la absolute bound of 13.3 x 10 
GeV. The relevant polar angles arc 6 ~ 99.0°, (f) ~ 82.6°, 
and the colatitude is x ^ -54.1°. while the neutrino energy 
is in the neighborhood of 20-60 MeV. Taking |p| = 40 
MeV as the representative energy, we estimate the max- 
imal sensitivity achieved to individual coefficients for 
Lorentz violation by setting all but one to zero and con- 
sidering separately any real and imaginary parts. The 
results of this calculation are displayed in Tables VIII 
and IX. 

Additional limits can be obtained from the recent 

study of Lorentz violation by the MiniBooNE collabo- 
ration [13]. This analysis placed bounds on amplitudes 
for both — > i/f, and 77^ — > Ve transitions with sidereal 
harmonics m = and m = 1. For fixed dimension d, the 
experiment gives 2a limits for v ^ u oscillations of 
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(4?)r+(4?)"^ 



< 
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GeV, 
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GeV. 



(101) 



For antineutrinos, the 2a limits are 

\pr'\{Am':'''+{Am-f < 



2.6 X 10" 


20 


GeV, 


3.7 X 10" 


20 


GeV, 


3.9 X 10" 


-20 


GeV. 



(102) 



The MiniBooNE analysis also assumes real Lorentz- 
violating terms in /loso which implies we must impose 
conditions on the coefficients as before. The polar angles 
for the beam direction are 6 ~ 89.8°, (/> ~ 180°, and the 
colatitude is x — 48.2°. The average neutrino energy is 
0.36 GeV, while the average antineutrino energy is 0.60 
GeV. Proceeding as above, we extract estimated maxi- 
mal sensitivities to individual coefficients for Lorentz vi- 
olation. The neutrino and antineutrino results obtained 
in this way are also compiled in Tables VIII and IX. 

The two tables contain many first limits on neutrino 
coefficients for Lorentz violation with d ~ 2 and with 
5 < d < 10. Several options exist for achieving improve- 
ments and extensions of these results. One possibility 
would be to reanalyze the LSND and MiniBooNE data 
for higher harmonics according to Eq. (98). This would 
generate first constraints on many additional coefficients. 
Another possibility is to use existing data from other 
short-baseline experiments. For example, data from the 
MINOS near detector have already been used to con- 
strain Lorentz violation within the short-baseline approx- 
imation [10], with sensitivities at the level of 10"^*^ GeV 



to two coefficients a 



for CPT-odd violation and at 



10 to seven coefficients 



for CPT-cven viola- 



tion. The same data could be analyzed for harmonics 
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TABLE VIII: Maximal attained sensitivities on the modulus of coefficients with d <6 from LSND antineutrinos at la [8] and 
from MiniBooNE neutrinos (MB) and antineutrinos (MB) at 2a [13]. The units are 10"^° GeV*"'^. 



to yield measurements of eoefficients with other values 
of d. Future short-baseline experiments such as the re- 
cent DAE5ALUS proposal [118] would also offer interest- 
ing possibilities for searching for Lorentz violation along 
these lines, as would analyses of data from reactor exper- 
iments such as Double Chooz [119], Daya Bay [120], and 
RENO [121]. 



B. Perturbative Lorentz violation 



The short-baseline approximation is appropriate for 
null experiments when little or no neutrino oscillation is 

detected. However, many experiments observe significant 
oscillations. In this general theoretical analysis is 

challenging. One approach is via model building, which 
in the context of Lorentz violation involves designing spe- 
cial Lorentz-violating models that can qualitatively re- 
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coefficient 


LSND 


MB 


MB 


0(6) 

= (8) 
.(10) 


2.9 
36 
910 

2.3 X 10'' 
5.7 X 10'' 

1.4 X lO'^ 
3.6 X 10^ 
8.9 X 10^ 
2.2 X 10" 


3 

4.2 

12 

32 

90 
250 
690 
1900 
5400 


3.1 
2.6 
4.3 
7.2 
12 
20 
33 
56 
93 



TABLE IX: Maximal attained sensitivities on the modulus of 
isotropie coefficients with d < 10 from LSND antineutrinos at 
la [8] and from MiniBooNE neutrinos (MB) and antineutri- 
nos (MB) at 2fT [13]. The units are 10"^° GeV"'^. 

produce the observed global features of oscillations, per- 
haps including also one or more of the neutrino anoma- 
lies, using only a few parameters [18-22]. Another strat- 
egy assumes oscillations are primarily due to the mass 
matrix and treats Lorcntz violation as a small pertur- 
bation, seeking to identify or constrain small deviations 
from the conventional picture that may indicate Lorentz 
and CPT violation. For the renormalizable limit of the 
SME, this approach is presented in Ref. [17]. 

In the present subsection, we gcineralizc the perturba- 
tive treatment to include Lorentz-violating operators of 
arbitrary dimension. Our primary focus is on beam ex- 
periments, in which the signal dependence on propaga- 
tion direction typically manifests as sidereal time depen- 
dence, but the basic analysis is applicable to other situa- 
tions such as the azimuthal signal dependence used in the 
recent search for Lorentz violation by the IceCube collab- 
oration [12]. For example, future searches with IceCube 
and Super-Kamiokande [122] could adopt the methods 
presented here to extract competitive limits on a variety 
of SME coefficients for Lorcntz and CPT violation. Other 
possible applications include searches for anomalous an- 
nual variations in solar-neutrino oscillations beyond those 
due to the Earth's orbital eccentricity. 

The perturbative approach is based on time-dependent 
perturbation theory, which expands the time-evolution 
operator S{t) in powers of the Lorentz-violating part Sh 
of the effective hamiltonian (51). This requires first spec- 
ifying the unperturbed system, which implies adopting 
values for the conventional mass-squared differences and 
mixing angles. To enable a straightforward match to pre- 
vious results, we adopt here the notation of Ref. [17] with 
upper-case indices AB indicating components of the full 
6x6 effective hamiltonian, and lower-case unbarred and 
barred indices ab, ah, ah, ah indicating components of 
the 3x3 blocks. As before, unprimed indices refer to the 
flavor basis while primed indices indicate the diagonal 
energy basis. In this notation, we write the conventional 
energy eigenvalues as Ea', with CPT invariance guar- 



anteeing the condition Ea' = Ea'- The 6x6 Lorentz- 
invariant effective hamiltonian (/ieff)oA_B of Eq. (38) is 
diagonalized by the 6x6 mixing matrix 

where U„,i, = U*,t is the familiar 3x3 neutrino mixing 

a o a'b ^ 

matrix. Note that the unperturbed system could be de- 
fined to include, for example, the effects of matter for 
neutrinos propagating through the Earth. These would 
alter the energies Ea' and the mixing matrix Ua'b- Since 
the interactions of neutrinos with matter are described by 

f3)i 

coefficients {a)^' )ab for CPT-odd Lorcntz violation [15], 
the oscillations of neutrinos and antineutrinos in matter 
involve different energy spectra and mixing matrices. 

The perturbation calculation generates an expansion 
of the oscillation probabilities in powers of the Lorentz- 
violating perturbation (44) with components SKab- The 
expansion is [17] 

where 

p(0) _ I 0(0) 1 2 o(0)_V^rr* J J -zE^,t 

(105) 

are the unperturbed oscillation probability and time- 
evolution operator. The flrst- and second-order pertur- 
bations take the form 

p!^^.a = ^ti^[sf;n%), 

Pi^U = -i^Re(5i°r^i^i)+t^|wW|^ (106) 
with 

^is = Y.^M%)cD6hcD, 

CD 

^aI = J2 {M'^AB)cDEF6hcD6hEF, (107) 
CDEF 

where {Mab)cd and {Mab)cdef are factors depend- 
ing only on unperturbed quantities and the experimen- 
tal setup. They are given explicitly by Eqs. (32) and 
(34) of Ref. [17]. Note that the first-order probabil- 
ity pI^b^va vanishes when the conventional transition 
is zero, = 0. Under these conditions, the second- 

(2) 

order probability Pi)J^VA governs the dominant Lorentz- 
violating effects. Examples of this arise when the only 
nonzero coefficients are of Majorana type and hence lie 
in the off-diagonal blocks of J/i, causing mixing between 
neutrinos and antineutrinos. 

1. First-order perturbation 

The first-order perturbations are governed by 'H^^. 

Since the unperturbed system contains no neutrino- 
antineutrino mixing, only "H^^' and its CPT conjugate 
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■^ij^ contribute to the first-order probabilities. More- 
over, only 6hab enters the expression for hI^^ because 
Ua'a is block diagonal. Analogous results hold for the 
other blocks of "H^^ . Note that all four blocks appear in 
the second-order probabilities. 



We can construct explicit expressions for the four 
blocks of 'Hab ill terms of spherical coefficients. For ex- 
ample, for the neutrino- neutrino block we obtain 



ce 
djm 



(108) 



To analyze a given beam experiment, it is convenient to 
adopt as before the standard laboratory frame with the x 
axis pointing south, y pointing east, and z pointing verti- 
cally upwards [112]. In this frame, the beam direction is 
constant and the coefficients for Lorentz violation acquire 
time dependence due to the rotation of the Earth. To ex- 
press the result in terms of coefficients in the canonical 
Sun-centered frame, we adopt the strategy employed in 
the short-baseline case in the previous subsection. The 
factors (A^Q6'')ce and the beam direction p^^^ are deter- 
mined in the laboratory frame. The time dependence is 
revealed by using the Wigner matrices to rotate the coef- 
ficients for Lorentz violation to the Sun-centered frame. 

Implementing this procedure for all four blocks of H^ab 
gives 



ay(l) 


— |p|d-3gimw®Tg 

dva 


((4'')):^-(4^)):^), 


W« 

ab 


dm 




ab 


— |p|d-3gimw®T(j 
dm 




^ab 


— |p|d-3gimw®Tg 
dm 





The various amplitudes appearing in these equations take 



the compact forms 



mm 

mi 
ml 
mt 
mt 
ml 
mz 



cej 



(<^)^^ce 



{d)\ce 
eff / jm' 



cej 



^oAf,m{Ml^^h,{ai'^);i, 



E oAfjm {M 

cej 



b >' 



(1)^ 



cej 

cej 



-e (H, 



eff , 



jm' 



T.-^^^rn{Mg%{-ir[{giX-m)] 



cej 

E 

cej 



[\^eS Jj(-m)i 



(110) 



In these expressions, the factors s-A/jm are defined by Eq. 
(93), as before. 

The above results can be directly applied to searches 
for Lorentz violation with beam experiments. For exam- 
ple, consider a study of Lorentz violation with — >■ or 
t^n T^e oscillations in a long-baseline experiment. The 
explicit numerical values of the relevant factors (A^i][?)a6 

and (A1g]x )a6 for a variety of long baseline experiments 
are given in Table I of Ref. [17]. For any given experi- 
ment, the beam polar angles and the colatitude can be 
used to determine the direction factors sl^jm via Eq. (93). 
Substituting the results into Eq. (110) gives the ampli- 

tudes {^r:- mZ {^)Z -levant 

for the first order probabilities in terms of the effective 



(d)\ab 



^eff , 



jm 



and 



{d)\ab 



^eff . 



) . . For a cho- 

/ jm 



spherical coefficients 

sen mass dimension d and a specific sidereal harmonic m, 
the contribution to the first-order probability for v^^ — >■ 
mixing is then given by 



sifimz-sirm 



+lm 



efj.-\ 



efj, 



sin mco^T^ 
cos mto^T^^ . 

(Ill) 



The corresponding result for i^^ — i>e mixing is immedi- 
ately obtained by changing the sign of the CPT-odd am- 
plitude (^A^'a ^y^ and replacing with ep, in the above 
expression. 

The result (111) shows that the given long-baseline ex- 
periment can use measurements of sidereal variations at 

various harmonics m to place constraints on linear combi- 
nations of the real and imaginary parts of the coefficients 
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(ai'S)'^.^ and (cii})°'.^ for Lorcntz violation. If a definite 
signal is found, a variety of independent experimental 
configurations may be required to identify which individ- 
ual coefficients are nonzero. In the absence of a signal, 
a useful approach is to take one coefficient nonzero at a 
time. The experimental bound can then be interpreted as 
a maximal attained sensitivity to each coefficient in turn, 
yielding tables analogous to Tables VIII and IX obtained 
above for the short-baseline approximation. Reporting 
the results in this way facilitates the direct comparison 
of different experiments and provides an effective guide 
to constraints on theoretical model building [3] . 



mixing channel can be expected to yield clean constraints 
on the coefiicients (.g|,fl')"^ and 

In principle, another possibility for achieving sensitiv- 
ity to the coefficients (ffoff )jm ^^'^ (^cff )jm *° 
sider instead second-order effects in neutrino-neutrino 
and antineutrino-antineutrino oscillations. The corre- 
sponding probabilities 



p(2) 



(1)|2 



pj^'l^ = -t'Re(s^fn^'A+t'\n['M' (iis) 

i^b^^a y ab ab J I ao I ^ 



2. Second-order perturbation 

The second-order probability becomes important when 
conventional oscillations arc negligible. One example is 
mixing between neutrinos and antineutrinos. Also, for 
certain energy and baseline combinations, v ^ u and 
77 — > 17 transitions may involve suppressed first-order ef- 
fects. 

For neutrino-antineutrino mixing, the first-order prob- 
abilities vanish. The second-order probabilities take the 
simple form 



p(2) _ ^ 



Expressions for H^^^ and T-L^^^ 



(1)1 

ab I 



pp) 



= t"H 



ab 



\ (112) 



1) 



are given in Eq. (109), 
with the corresponding amplitudes expanded in terms of 

the effective spherical coefficients (ffeff ) ^^'^ (^es)'jm 
in Eq. (110). These results can be applied to searches 
for Lorentz violation in neutrino-antineutrino mixing, us- 
ing a procedure similar to that outlined in the previ- 
ous subsection for long-baseline experiments on neutrino- 
neutrino and antineutrino-antineutrino mixings. Since 
neutrino-antineutrino oscillations cannot arise from the 
coefficients [o-'^)'^^ and (c^g )"^, direct studies of this 

I 



depend on the second-order combinations 



n 



(2). 
ab 



n 



(2)^ 
ab 



Yl [i^ab^cefgSKeShfg 
cefg 



SgfgShceShfg 



i^ni)BefgShceShf-g 



cefg 



(114) 

This line of attack is more complicated because the com- 
binations (114) arc quadratic in coefficients for Lorentz 
violation and because they contain all four 3x3 blocks 
of 6h. Although the diagonal blocks Shab and ^/i^g also 
contribute to the first-order probabilities, which implies 
their appearance in Eq. (114) represents an effect only 
at subleading order, the combinations of coefficients ap- 
pearing in the probabilities may be distinct. As a conse- 
quence, the contributions of the diagonal blocks cannot 
be omitted in an exact treatment. 

Nonetheless, to gain intuition about the leading-order 
contributions to Eq. (114) arising from (.9^ff ) ; and 



id)) 
eS ) 



jm' 



we can choose to focus on the off-diagonal 



blocks 5/1^5, Shab while disregarding the contributions 
from Shab, Sh^i- For we thereby obtain 
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(115) 



cefg jj' 
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These equations explicitly confirni that data from 
neutrino-neutrino mixing also contain information about 

the coefEcients (g^g-*)"^ and (iJ^g')"^^. An exact 
treatment would generate constraints on complicated 
quadratic combinations of all four sets of coefficients 

(Jd)Yib ( {d)\ab (TT{d)\ab 

In the absence of a definitive signal in the data, an 
interesting approach is to take only one coefficient for 
Lorentz violation to be nonzero at a time. Sidereal vari- 
ations can be used to determine maximal sensitivities 
to individual cartesian coefficients for Lorentz violation 
[17]. In the present context, considering a single effective 
spherical coefficient implies working with fixed values of 
d, j, m. The only terms that survive in the above ex- 
pressions then have a unique value of m = — m', and so 
no sidereal variations appear. Moreover, the cross terms 



of TifJ therefore takes the simple form 



between (0!^ and (H^f))'^.^ vanish. The exact form 



n 



ab 



S9 Jm{-m)' 



d odd, 
d even. 



(116) 



Note that the contributions involving odd d come from 

operators that are CPT even, while those involving even 
d come from operators that are CPT odd. Also, only 
even powcirs of the energy appear. The unconventional 
energy dependence offers a potential experimental handle 
for studies of (9^3)°^^ and (H^g)'^^. Another option 
arises for null experiments producing an absolute bound 
on the total oscillation, which can use the above simple 
expression to extract maximal attained sensitivities to 
individual components (,9cff and (^^og) 



C. Two-flavor maximal mixing 

The special case of two-flavor mixing with the max- 
imal mixing angle of 45° offers a simple example with 
phenomenological relevance. In the Si^SM, the smaller of 
the two mass-squared differences becomes less important 
at higher energies, so oscillations become dominated by 
the larger mass-squared difference. This situation applies 
to the energy regime relevant for most experiments with 
atmospheric and accelerator neutrinos. Observations in- 
dicate the 3z^SM angle 9is is small and 6*23 — 45°, which 
implies oscillations between i/^ and can be well ap- 
proximated by the two-flavor maximal-mixing scenario. 

In this limit, the only nonzero flrst-order transition 
probabilities are P^P-^^^ = Pi^L).„ , given by 



p!)IUu^ = L sin(Am2L/2E)Re((5V) 
= sin(ATO^L/2£') 

(d)\{t^T) 



djm 



where the symmetrized combinations are deflned by 

{"'isY-IL '' = ("-cfi- + ("'off 

and obey the conjugation relations 
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Ad}\if^'r) 
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m(id)\itiT) 



ji-m)' 



(119) 



(1) 



The survival probabilities are given by Pv^ 



1 - 

r t^u, 



(1) 



For antineutrinos, the transi- 



We remark in passing that other types of experi- 
ments may also yield direct sensitivity to (ffeff )j,„ arid Pj^^Lj 

[h'^^Y^^^^. For example, neutrinoless double-beta de- tion probability P^^^'^-p and the survival probabihties 
cay can arise in the presence of Majorana-like couplings 
that mix neutrinos and antineutrinos, and so studies of 
this process could provide competitive measurements for 
Lorentz violation mediated by Majorana operators. A 
treatment of these and related possibilities would be of 
deflnite interest but lies outside our present scope. 



= P^^^-^ are obtained by changing the sign of 

the coefficient {a^3)^^'^\ 

In terms of coefScients in the Sun-centered frame, the 
transition probability acquires a sidereal time depen- 
dence, taking the form 



Pl%.. = |isin(Am2L/2E)^|p|<^-3[oW;o(ai^ff^);r^-o^,o(cS)) 
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(120) 
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where the factors oA/^m arc given by Eq. (93), as be- 
fore. The combinations of coefficients for Lorcntz viola- 
tion that can be measured in this scenario arc therefore 
the real and imaginary parts of o-^m(aeff ) -^-^ 



^ J"" \ en /jm 

This scenario also provides a simple framework for 
studying the effects of CPT violation. For example, con- 
sider the CPT asymmetry 



(121) 



Assuming identical energies, baselines, and beam direc- 
tions, this asymmetry vanishes when CPT is conserved. 
A nonzero experimental measurement of A'^^^ would 
therefore provide evidence of CPT violation. We remark 
in passing that a zero measurement would fail to prove 
CPT symmetry because CPT-violating models exist for 
which A^j^'^ vanishes [15]. 

For maximal two-flavor mixing, the asymmetry A^^"^ 
can be expressed compactly in terms of coefficients for 
Lorentz violation in the Sun-centered frame. Assum- 
ing the conventional zeroth-order transition probability 
p!)%^^ = sir?{I^m?L/AE) is large compared to CPT- 
violating effects, appearance experiments are appropriate 
and the CPT asymmetry takes the form 



aCPT 



L cot{Am^ L/AE) 

djm 



-3 ^imujQ 



(122) 



If instead the survival probability Pi^Ui^^ is large com- 
pared to CPT-violating effects, then disappearance ex- 
periments are useful and the relevant CPT asymmetry 
is 



A 



CPT 



djm 



f 3m 



(123) 



The above two asymmetries contain the same coeffi- 
cients for Lorentz violation but apply in complementary 
regimes, so at least one is applicable to any experiment. 
They incorporate direction-dependent effects and varia- 
tions with sidereal time as a consequence of the Lorentz 
violation that accompanies CPT violation. Similar ef- 
fects accompany CPT violation in the oscillations of neu- 
tral mesons [123], and indeed the above expressions are 
closely related to the corresponding CPT asymmetries 
for mesons [124]. 

The CPT asymmetries contain coefficients for Lorentz 
and CPT violation with m = that produce effects inde- 
pendent of sidereal time. One simple way to extract these 
coefficients, already used in the meson context [124], is 
to measure the time-averaged asymmetries A^^'^ and 
"^uif"^- -'^O'^ given energy, baseline, and beam direction. 



these asymmetries take the form 

« Lcoi{Am'L/AE)Y,\p\'-\MA<^^^)t\ 



dj 



AC^^ « -Ltan{Am'^L/AE)Y,\pr\^f,,{a^^f'^^ 



dj 



jo ■ 

(124) 



The presence of the factors o-A/jo shows that these ex- 
pressions depend on the beam direction despite the time 
averaging. As a result, distinct experiments can be ex- 
pected to have different sensitivities to different coeffi- 
cients for Lorentz and CPT violation. Over the range of 
existing and planned long-baseline experiments, a given 
factor oA/jo can change sign and can vary by more than 
an order of magnitude, so the difference in attained sen- 
sitivities can be substantial. 



VII. APPLICATIONS TO KINEMATICS 

Neutrino oscillations can yield observable signals of 
Lorentz violation because they involve comparing the 
propagation of one neutrino flavor against another. An- 
other possibility for detecting physical effects of Lorentz 
violation is to compare neutrino propagation to the prop- 
agation of a different kind of particle. Kinematic tests 
of this kind come in several varieties. One conceptually 
straightforward approach is to measure the difference in 
the times of flight of neutrinos and photons or other par- 
ticles. A more subtle possibility is to study decay pro- 
cesses involving neutrinos, which can be modified when 
the dispersion relations of neutrinos and other species 
differ in their Lorentz properties. The effects on decays 
can be striking, with certain processes becoming forbid- 
den or allowed according to the energies of the particles 
involved. 

In this section, several kinds of kinematic tests are con- 
sidered. To focus the discussion, we assume oscillations 
are negligible or zero. This implies working within the 
context of the oscillation- free models described in Sees. 
V C 2 and V E 3. Both generic oscillation-free models and 
the isotropic flavor-blind limit are treated. We consider 
time-of-flight measurements, threshold effects in pion and 
kaon decays, and Cerenkov radiation, using existing data 
to obtain explicit results for coefficients for Lorentz vio- 
lation. 



A. Time-of-flight measurements 

Time-of-flight experiments compare the group veloc- 
ity of neutrinos with that of photons or other particles. 
Here, we work with four measurements involving time-of- 
flight comparisons with photons: the recent OPERA re- 
sult [14], the prior MINOS bound [125], early constraints 
from experiments at Fermilab [126], and limits from the 
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TABLE X: Numerical values of direction factors oMjm for some long-baseline experiments. Values are given to two decimal 
places. An explicit 0.00 indicates a small nonzero value. 
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9.1 ±5.4 X 10"® 


< 


4.6 


X 


10" 


8 










40 


24.2 ± 3.3t2;6 X 10"* 


5.2 ±3.0 X 10"® 
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TABLE XI: Single-coefficient measurements and modulus bounds from accelerator time-of-flight experiments. Units are 
GeV""". 
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supernova SN1987A [127]. Lorentz violation in the pho- 
ton sector involving operators of both renormalizable and 
nonrenormalizable dimensions is tightly constrained [3]. 
For definiteness, we assume below a conventional photon 
dispersion relation, with any Lorentz violation confined 
to the neutrino sector. 



1. Generic case 



„of 



In generic oscillation- free models, the group velocity 

for neutrinos propagating in direction p is given by 
Eq. (74). For astrophysical neutrinos, this result is di- 
rectly applicable. However, for neutrinos in a beam ex- 
periment, sidereal variations in the signal are induced by 
the rotation of the Earth. These can conveniently be 
handled via the methods discussed in Sec. VIA. Work- 
ing as before in the standard laboratory frame [112], the 
Wigner rotation matrices can be used to display the side- 
real variations explicitly and to express the; group veloc- 
ity in terms of spherical coeflScients for Lorentz violation 
in the canonical Sun-centered frame [111]. The effect of 
this procedure on the group velocity (74) is to perform 
the substitution Yjm{p) -> e*""®^®o-A/'jTO, where the fac- 
tor o-A/jm is defined in Eq. (93). We thereby find for beam 
experiments the neutrino group velocity 



= 1 - 



djm 



x[(-i?),>„-(ei?U>(125) 



which displays explicitly the dependence on the sidereal 
rotation frequency and the local sidereal time T^. 
We remind the reader that for antineutrinos the sign of 
the coefficients {a'^^^)j^ changes in all expressions for the 
group velocity. 

The OPERA collaboration reported a difference be- 
tween the speed of light and the speed of muon neutrinos 
of Sv = 2.37 ± 0.32;°;34 ^ ji4j^ Averaging over 

sidereal time, this velocity defect yields the condition 



dj 



Jd)) 



Jd)) 



of JjO V-of JjO 

= 2.37 ± 0.321^:^^ X 10"^ (126) 

on the oscillation-free spherical coefficients for Lorentz 
violation. Here, oMjo is the directional factor for the 
OPERA beam, which must be computed using the beam 
angles = 86.7°, = 52.4° at the colatitude x = 47.5° 
of the detector. Numerical values for relevant values of 
oJ\fjo for OPERA are listed in Table X. 

To gain intuition and for purposes of comparison with 
other experiments, we can extract from the condition 
(126) a set of constraints on individual oscillation-free 
coeflacients for Lorentz violation under the assumption 
that only one is nonzero at a time. Using for |p| the av- 
erage energy of the beam, ~ 17 GeV, the resulting 
constraints are listed in Table XI for dimensions d < 6. 



The shape of the observed neutrino spectrum matches 
the expected form to a high degree, which implies little or 
no dispersion in the group velocity at OPERA energies. 
This is reflected in the comparison of neutrino group ve- 
locities at low ((IpI) ~ 13.9 GeV) and high ((|p|) ~ 42.9 
GeV) energies. The arrival-time difference reported by 
OPERA is 14.0 ±26.2 ns, which translates into a velocity 
difference A?; between the two datasets of approximately 
Av ~ 6 ± 11 X 10~^. As an illustration, we can use this 
to place a comparatively reliable constraint on dimension 
d = 5 operators for Lorentz and CPT violation, assum- 
ing Lorentz-violating operators at other values of d are 
negligible. We thereby find the condition 



J^o^jo (aifO.o - 10 ± 19 X 10-8 GeV- 



(127) 



on a beam-dependent combination of coefficients with 
(i = 5. The lack of dispersion could in principle also 
be used to bound coefficients for d > 5, but obtaining 
reliable constraints requires access to more detailed in- 
formation about the observed energy spectrum. 

A prior time-of-flight experiment by the MINOS col- 
laboration also measured the group velocity for the muon 
neutrino compared to the speed of light, with the result 
6v = 5.1 ± 2.9 X 10~^ [125]. For generic oscillation-free 
coefficients, we obtain the condition 



J:{d-3)\pf-\^^,o[{a'^^).^-{cif)^^ 

dj 



5.1 ±2.9 X 10~ 



(128) 



where we have again averaged over sidereal time. The 
numerical values of the directional factors oMjo for MI- 
NOS arc listed in Table X. They are computed using 
the polar angles = 86.7°, (j) = 203.9° of the beam and 
the colatitude x = 42.2° of the detector. Adopting the 
average beam energy as (|p|) = 3 GeV, we can extract 
single-coefficient constraints taken one at a time. For 
dimensions d< 6, the results are listed in Table XI. 

An older experiment at Fermilab [126] reported bounds 
of \6v\ < 4x 10^^ at 95% C.L. using both muon neutrinos 
and muon antineutrinos. Averaging over sidereal time, 
this gives the two conditions 



)jO 



^(rf-3)|p|''-\^,o[±(al?).3-(ci?)^ 

dj 

< 4 X 10"-\ (129) 



The experiment also yielded a limit on the difference be- 
tween the neutrino and antineutrino group velocities of 
7 X 10~^. In the present context, this generates the limit 

|$^(rf-3)|pr%AA,-o(ai?).„| <3.5xl0-^ (130) 

dj 

on a combination of coefficients for GPT-odd Lorentz vi- 
olation. As before, we can place constraints on individ- 
ual coefficients taken one at a time. We adopt the esti- 
mated beam angles ~ 90°, </> ~ 140° and the colatitude 
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TABLE XII: Single-coefficient modulus bounds from time-of-ffight and dispersion of SN1987A antineutrinos. Units are GeV* 



X = 48.2° for the detector. In these experiments the en- 
ergies ranged from 30 to 200 GeV, so for our calculations 
we take the conservative value \p\ = 30 GeV. The results 
for dimensions d < 6 are listed in Table XI. 

Observations of the antineutrino burst from supernova 
SN1987A lead to a conservative bound on the difference 
between the speed of light and the speed of antineutrinos 
of \Sv\ < 2 X 10^^ [127]. The large propagation distance 
implies the electron antineutrinos produced at the source 
oscillated many times during their trip to the Earth. No 
sidereal effects occur in this case, so we obtain the con- 
dition 



djm 



{d~'6)\pf-'Y,m (al?),„ + (cin 



< 2 X 10" 



(131) 

involving both isotropic and anisotropic coefficients. This 

means we can also extract here individual constraints on 
coefficients with nonzero values of m. In the Sun-centered 
frame, the propagation direction is given by the polar 
angles 9 = 20.7°, (j) = 263.9°. The observed energies 
range between about 7.5 and 40 MeV, so we adopt the 



conservative value \p\ = 10 MeV. The resulting single- 
coefficient time-of-flight constraints for values d < 6 are 
given in Table XII. 

The observed antineutrinos from SN1987A have a 
spread of energies. However, they all arrived within 
a time interval of about 10 s after traveling for about 
5 X 10^^ s, implying a maximum difference in speed 
5v < 2 X 10~^^ across the observed energies. This re- 
stricts the possible antineutrino dispersion and implies 
additional constraints on coefficients for Lorentz viola- 
tion, independent of the speed of light [15]. For an energy 
spread ranging from \pi\ to \p2\, the dispersion condition 
is 



djm 



{d-3)A{\pf-^)Y, 



|d-4 



jm 



jm 



< 2 X 10"^^ (132) 

Using the same 



where Adpl''-^) = \p2r— - |Pi 
values of and </> as before and adopting the conservative 
choices \pi\ = 10 MeV and jpsj = 20 MeV yields the 
single-coefficient dispersion constraints for values d < 6 
listed in Table XII. 
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15 



Fermilab 



< 5.8 X lO"'^ 

< 3.2 X 10 

< 2.4 X 10 



1-10 



1-13 



SN1987A 

< 2.0 X 10"^ 

< 1.0 X 10"^ 

< 6.7 X 10"® 

< 5.0 X 10"^ 

< 4.0 X 10"^ 

< 3.3 X 10° 

< 2.9 X 10^ 



SN1987A 



< 1.0 X 10"^° 

< 2.2 X 10"® 

< 7.1 X 10"® 

< 2.7 X 10"® 

< 1.1 X 10"* 

< 4.5 X 10"^ 



TABLE XIII: Time-of-fliglit and dispersion measurements and modulus bounds on isotropic oscillation-free coefficients. Units 

are GeV*"'*. 



In principle, direction-dependent constraints can also 
be extracted from beam experiments on the Earth [15, 
123]. The sidereal variations due to the Earth's rotation 
cause the neutrino group velocity to change with time. 
Note that the observed group velocity can oscillate be- 
tween superluminal and subhiminal values. The sidereal 
variations are encoded in the factor cxp{imLu^T:^) ap- 
pearing in the group velocity (125). For fixed d, we can 
introduce complex amplitudes associated with a particu- 
lar sidereal harmonic as 

= E(^-3)0A/;>«(ai?),„, 

j 

i^i'^)m = E(^-3)oAA,™(ci?).„, (133) 

j 

which obey (^if])^ = (-4ifl)_^. The neutrino velocity 

defect in oscillation-free models can then be written 
as the expression 



(^i^'0„-(4^'')J. (134) 



dm 



The antineutrino velocity defect takes the same form 

except for the CPT-induced sign change of (a'^^)^. 

The above sidereal variations of the group velocity 
can be used to extract constraints from existing and fu- 
ture time-of-flight data obtained in beam experiments. 
The beam direction relative to the Earth fixes the oAfjm 
factors and hence determines the relevant linear com- 
binations of coefficients for Lorentz violation. To illus- 
trate the variations in these combinations as a function 
of beam direction, we list in Table X the oA/jm factors 
for the long-baseline experiments MINOS, OPERA, and 
T2K [128]. For the MINOS and OPERA experiments, 
we use the polar angles and detector colatitudes given 
above. For the T2K experiment, we adopt the beam 
angles 6 ~ 88.7°, (f) — 270° and take the detector co- 
latitude as X = 53.6°. Note that the reality condition 
(94) implies oN*^ = (— l)"'oA/^(-m) j which can be used 
to compute the factors o^jm for negative m values. The 
values in Table X can be inserted in the expressions (133) 
and thereby into the velocity defect (134) to extract lim- 
its on the oscillation-free coefficients in several time-of- 
flight experiments. Relevant searches include ones with 



MINOS, OPERA, and T2K, as well as Borexino [129] 
and ICARUS [130], both of which have the same factors 
oMjm as OPERA. 



2. Isotropic case 

The discussion in the preceding subsection holds for 
generic oscillation-free models. In contrast, for isotropic 
oscillation-free models, the physics is independent of the 
beam direction and so no sidereal variations arise. The 
neutrino group velocity v for isotropic oscillation-free 
models is instead given by Eq. (87), which applies equally 
to astrophysical neutrinos and neutrinos in beam exper- 
iments. As before, the antineutrino group velocity v dif- 
fers by a change of sign for the coefficients o/'^\ 

In the isotropic case, the OPERA measurement [14] 
yields 

J2{d - 3)|p|''"*(a('^) - c^'^^) = 2.37 ± 0.3215]^^ x 10"^ 

d 

(135) 

The lack of dispersion in the neutrino pulse yields a com- 
paratively reliable constraint on dimension-5 operators 
for isotropic Lorentz and CPT violation of 



^(5) 



10 ± 19 X 10"* GeV" 



(136) 



Using the MINOS result [125], we obtain the condition 
"^{d-3)\p\'^-^{a^'^^ -c^'^^) = 5.1±2.9xl0"5 (137) 

d 

on isotropic oscillation-free coefficients. Analogously, the 
older Fermilab results [126] yield the bound 

I ^{d - 3)|p|'^"4(±a('*) - c^'^^) < 4 X 10"^ (138) 



along with the constraint 

l^id-SM'-'a^"^ 



< 3.5 X 10"^ 



(139) 



on CPT-odd effects. Finally, the SN1987A observations 
[127] give the time-of-flight bound 

\"^{d - 3)\p\'^-^ {a^'^^ + c^'^^) <2xl0"9 (140) 



28 



and the dispersion bound 

|^(d-3)A(|p|''-^)(a('^)+c('^); 



< 2 X 10"^^ (141) 



Under the assumption that only one term is nonzero 
at a time, we can use the above conditions to extract 
Umits on isotropic oscillation-free coefficients for Lorentz 
violation. For Lorentz-violating operators of dimension 
d < 10, the results are displayed in Table Xlll. Future 
long-baseline experiments and astrophysical observations 
can be expected to improve these constraints. 



B. Threshold efFects 

The decay processes 7r+ — >■ ^+ + 1'^ and — >■ /x+ -|- 
are the dominant sources of muon neutrinos in most ex- 
periments. In the presence of unconventional dispersion 
relations for neutrinos, these decays can exhibit strik- 
ing threshold effects [51, 70-72, 80-82, 96, 99, 131]. In 
this subsection, we use threshold effects in pion and kaon 
decays to obtain additional constraints on the spherical 
coefficients under the assumption that any Lorentz vio- 
lation appears only in the neutrino sector. 

For pion decay, let k be the pion momentum and p be 
the neutrino momentum. Energy-momentum conserva- 
tion implies the neutrino energy E{p) obeys 

E{p) = ^Ml+k^-^M^ + ik-pf 



< ^M2+fe2_^M2 + (|fe|-|p|)2 

< V'AM2Tp2, (142) 

where AM = — M^. The first inequality is obtained 
by taking p parallel to k to maximize the allowed en- 
ergy. The last relation is a reverse triangle inequality. 
Similar relations hold for kaon decay, with the replace- 
ment Af^ Mk- 

If the neutrino energy E{p) is Lorentz invariant, the 
above inequalities are satisfied. However, a Lorentz- 
violating neutrino dispersion relation can cause the in- 
equalities to fail above some threshold energy, in which 
case the decay becomes forbidden by energy-momentum 
conservation. Experimental observations of high-energy 
muon neutrinos therefore constrain modified dispersion 
relations. To identify the conditions resulting from the 
energy threshold, we can write 



E{p) = Eo{p)+5E{p), 



(143) 



where Eo{p) is the conventional neutrino energy and 
SE{p) is the Lorentz-violating contribution. The reverse 
triangle inequality then gives 



SE{p) < 



AM2 



2En 



< 



AM2 



(144) 



The latter inequality can be applied to obtain explicit 
one-sided constraints on Lorentz violation. 

For the generic oscillation- free coefficients, the expres- 
sion (73) for the neutrino energy E°^ and the CPT con- 
jugate for E^^ yield 

= Y: \p\'-' YMp)1± - (4?), J . (145) 

djm 

We therefore obtain two one-sided constraints. 



E 

djm 



p 



d-2 



Yjm{p) 



< iAM^ 



(146) 

Notice that these bounds depend on the direction of the 
neutrino or antineutrino propagation. 

The IceCube collaboration observes atmospheric neu- 
trinos at high energies up to about 400 TeV [132]. Using 
the inequality (146), data from this experiment could 
be used to search for directional effects involving a re- 
duced muon-neutrino flux depending on the polar angles 
{9,(f>). This study would be qualitatively different from 
the search for Lorentz violation in oscillations of atmo- 
spheric neutrinos already published by the IceCube col- 
laboration [12]. It would provide sensitivity to distinct 
coefficients for Lorentz violation. 

For the isotropic oscillation- free case, the neutrino en- 
ergy Ev given in Eq. (86) and the corresponding antineu- 
trino energy E17 yield instead the two conditions 

SE = J2 \pf~^i±ci^'^^ - c^'^^), (147) 



which lead to the two one-sided bounds 



^|p|''-2(±a('')-c('')) < iAM^. 



(148) 



For pion decays, the numerical value of the right-hand 
side is ^AM^ = 5.7 x lO"'^ GeV^ For kaon decays, 
the result is weaker by roughly two orders of magnitude, 
iAM|. = 7.5 X 10-2 GeV^. The kaon mode typically 
dominates at higher energies, but heavier mesons may 
contribute as well. 

Using the inequality (148) reveals that the IceCube 
observation of neutrinos up to about 400 TeV suffices by 
itself to place fairly robust limits on isotropic oscillation- 
free coefficients for Lorentz violation. We obtain 

^(400TeV)'^-2(±a('')-c('')) < 7.5x10-2 GeV^. (149) 

Since no significant deviation is seen in the total neutrino 
and antineutrino flux, we can reasonably assume this 
constraint applies independently to both muon neutrinos 

and antineutrinos. Taking only one spherical coefficient 
to be nonzero at a time therefore produces two kinds 
of constraints. A lower (negative) bound is obtained on 
coefficients for CPT-even Lorentz violation, while a two- 
sided bound emerges for the CPT-odd case. Table XIV 
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a(^'| 


< 7.3 X 10" 


30 


g(8) 
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-1.8 X 10" 


-35 






< 4.6 X 10" 


41 


^(10) 


> 


-1.1 X 10" 


-46 



TABLE XIV: Estimated bounds on isotropic oscillation-frcc 
coefficients from a thresliold analysis of IceCube data. Units 
are GeV*"''. 



lists the resulting estimated bounds on single isotropic 
oscillation-free coefficients for values d < 10. 

For the specific isotropic oscillation-free model with 
only the coefficient c^^^ nonzero, two bounds from thresh- 
old effects in meson decay have recently been obtained 
using the IceCube observation of neutrinos with energies 
up to 400 TeV. The first translates into the one-sided 
bound c^*' > -4 X 10~^^ [81], and the second into the 
one-sided bound c^^) > -lO^^^ [32]. The first bound is 
two orders of magnitude tighter than our result in Ta- 
ble XIV because it assumes IceCube neutrinos originate 
from pion decay. The second result is consistent with our 
analysis. All the results are many orders of magnitude 
more stringent than the value of c^^^ obtained from the 
OPERA result and Usted in the first row of Table XIII. 



C. Cerenkov radiation 

Another kinematic feature of Lorentz violation is the 
possibility that a neutrino undergoes Cerenkov radiation, 
emitting one or more particles [23, 63, 68, 78, 79, 84, 
86, 92, 93, 97-99, 133]. This can occur when the maxi- 
mum attainable velocity of the neutrino exceeds that of 
the emitted particles. In this subsection, we comment 
on some implications of neutrino Cerenkov radiation in 
the context of the present work, with foc;us on thresh- 
old effects and spectral distortions. Except where stated 
otherwise, Lorentz violation is assumed to be confined to 
the neutrino sector. 



1. Threshold effects 

Suppose a neutrino with high energy E(j)) and with an 
unconventional dispersion relation experiences Cerenkov- 
like emission to one or more particles. Potentially signifi- 
cant Cerenkov processes for neutrino energies up to some 
tens of GeV include neutrino splitting -^■ i^,, + v,. + T'e, 
electron-positron pair production j/^ — >• + + e~ , 
and photon decay — >■ t'p + 7. However, as the neu- 
trino energy increases other processes become important, 
such as muon pair production 1^1x^1^^+ /x"*" -|- tau 
pair production z/^ ^ i/^ -|- r"*" -|- t~ , emission of various 
hadron combinations, and ultimately even emission 
i^n ^ ^fi + Z'^ and Higgs emission i^fx ^ 1^^ + (j). 



Energy-momentum conservation for the radiation of a 
single particle of mass M and momentum k implies 

AE{p, k) = E(p) - E{p -k) = -|- k^. (150) 

Similarly, for emission of two particles of masses Mi , M2 
and momenta fci, k2, we can write 

AE{p,k) = E {p)-E{p -k) 

= ^JMi + ki + ^M| + kl 

> VM2 + fe2 = A£th(fe), ^^5^^ 

where M = Mi + M2 and fe = fei -|- fe2. We see that 

AEth{k) represents a threshold on the emitted energy, 
obtained by treating the radiated particles as a single 
composite particle. This result generalizes to larger num- 
bers of conventional decay products. Note, however, that 
it relies on Lorentz violation being confined to the orig- 
inal particle. In processes of this type, natural scenar- 
ios exist in which modified dispersion laws for the emit- 
ted species exclude the presence of a finite threshold and 
hence forbid Cerenkov radiation. 

A given Cerenkov process is forbidden when the change 
AE{p, fe) in neutrino energy remains below the emitted- 
energy threshold A£'th(fe) for all values of k. This condi- 
tion can be visualized graphically by plotting AE{p, k) 
and A£'th(fe) versus |fc|. The threshold curve for emission 
of particles with nonzero total mass takes the form of a 
conventional mass hyperbola, while for purely massless 
emission it is a conventional light cone. The Cerenkov 
decay is allowed if the neutrino curve passes above the 
threshold curve. 

Suppose the momentum k transferred is small. The 
energy emitted is then AE{p, fe) « fe • v{p), where v{p) 
is the neutrino group velocity. This is maximized when 
fe lies along the direction of v{p). At the threshold, we 
then find k^v"^ = M^ + fe^, which implies 

fe2 = __ (threshold). (152) 

We emphasize that this result for the onset of Cerenkov 
radiation is independent of the specific form of the neu- 
trino dispersion relation. 

The recent OPERA measurement of the muon- 
neutrino velocity defect S\v\ = 2.37 ± 0.32^0 24 x 10"^ 
[14] yields the threshold value 

]fe| ~ 145M (153) 

for neutrino Cerenkov radiation into species of total rest 
mass M. For example, the threshold for electron-positron 
emission is |fe| ~ 150 MeV, while the threshold for muon 
pair production is |fc| ~ 31 GeV. The OPERA neutrinos 
have average energies ~ 17 GeV, with a spectrum 
exceeding 40 GeV. The result (153) therefore shows in a 
model-independent way that the neutrinos observed by 
OPERA are above threshold for Cerenkov emission. The 
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same argument also implies that the atmospheric neu- 
trinos exceeding 100 TeV observed by the IceCube col- 
laboration [132] lie well above the threshold of |A;| ~ 13 
TeV for Cerenkov emission. Note that the occurrence 
of more complicated processes such as neutrino splitting 
can only serve to strengthen these conclusions. 



Z. Spectral distortion 

The neutrino energy loss to Cerenkov radiation im- 
plies a distortion in the energy spectrum of the neutrino 
beam that depends on the baseline. Nonobservation of 
this distortion therefore offers a potential basis for plac- 
ing constraints on Lorentz violation. One measure of 
relevance for this analysis is the energy loss per distance 
dE/dx, which can formally be obtained for arbitrary co- 
efBcients for Lorentz violation using existing techniques 
for evaluation of scattering and decay processes in the 
SME. A complete analysis requires determining contri- 
butions from all Cerenkov processes above threshold. 

For each Cerenkov process, care is required to account 
for two distinct kinds of Lorcntz-violating modifications 
to conventional results. One is kinematical effects in 
phase space, associated with the unconventional disper- 
sion relations. These are comparatively straightforward 
to treat. The other is changes to the matrix element, 
which arise in several ways and are more subtle to han- 
dle. In processes with external fermions like neutrino 
Cerenkov radiation, the basic spinor solutions can be 
unconventional because Lorentz violation generically im- 
plies that spin no longer commutes with the hamiltonian 
[6, 134]. Also, gauge invariance implies that Lorentz- 
violating neutrino dispersion relations come with uncon- 
ventional interactions, so vertices such as vvZ^ acquire 
modifications. Moreover, in realistic models at least part 
of each neutrino field lies in an electroweak doublet with 
a charged lepton, so Lorentz-violating neutrino proper- 
ties imply modifications to charged leptons. As a result, 
any Cerenkov process involving leptons such as electron- 
positron pair production acquires accompanying Lorentz- 
violating contributions. These may qualitatively change 
the physical behavior, in some circumstances even elimi- 
nating Cerenkov emission. 

As an example providing some insight, consider the 
case of electron-positron pair emission in the presence of 
neutrino-sector Lorentz violation. If the neutrino energy 
is suificicmtly above the threshold (153) for this process, 
which is true for OPERA energies and higher, the neu- 
trino and electron masses can be neglected. For simplic- 
ity, suppose the only Lorentz violation in neutrino prop- 
agation is of Dirac type. The contributions from neu- 
trino propagation are then controlled by the term —V^aft 
in the Wcyl hamiltonian (30), or equivalently by the ef- 
fective coefficients a^g and c^g in the hamiltonian (44). 
Including the effect of Lorentz violation in the interac- 
tion vertices and in the electron sector would require fur- 
ther development of the formalism presented in this work. 



which lies beyond our present scope. However, these ef- 
fects are of the same order in Lorentz violation as the 
eflfects on propagation, so neglecting them can plausibly 
be expected to yield results of the correct order of magni- 
tude except in special circumstances. As a consequence, 
although the assumption of Lorentz violation confined to 
the neutrino sector is strictly inappropriate for a com- 
plete and realistic analysis of Cerenkov radiation, we can 
nonetheless expect to obtain reasonable insight by eval- 
uating effects from modifications of the neutrino spinors 
and of the kinematics. 

Under these assumptions, the relevant contribution to 
the matrix element in unitary gauge is 



iM = 



-iy/2GFMl 



xu{k)-i^{2s'' - PL)v{k'), 



(154) 



where the incoming neutrino has momentum p, the out- 
going one has momentum p' , the electron and positron 
momenta are k and k' , and s = sin 9w- To ensure valid- 
ity of the results at high neutrino energies above the 
pole, we keep the Z'^ propagator instead of adopting the 
four-Fermi approximation. 

In determining the square of the matrix element, the 
sum over electron spins yields the conventional result by 
assumption. However, the neutrino spin sum is modified. 
The neutrino spinors obey the modified Wcyl equation 
q ■a(j)L = 0, where q'^ = — is an effective lightlike 
momentum satisfying 



qo = E-Vl = \q\ « \p\ - 



P-Vl 

\P\ ' 



q/qo = {l,q) « ihp - V l/\p\ + PP ■ V l /IpI""). (155) 

Adopting the usual normalization for the spinors implies 

(l>L<P^L — -^(1 ~ 9 ■ ^) = ^1 ■ <^/Qo- Using this result, we 
obtain 

V- I |2 ^ 32GlM^{l-4s^ + 8s'')q-kq'-k' f EpE^A 
' [(fc + fcO^-M|P Uo9W' 

(156) 

where the factor in parentheses results from the modified 
neutrino spinors. 

The energy loss per distance is given by the integral 



dE _ _ f 



d^k 



d^k' 



f2Ep, {2TTf2Ek {2w)^2Ek' 
(^fc + Ek'){27r)'6Hp -p'-k'-k)j:\M\^ 
2Ep 

— J sin e sin 9' ded(f)de' d(t>' d\K\ 

k^k^k!'^ d\t^\ q-kq'-k' 
(k2 - M|)2 OkO qoEkq'oEk' 



(157) 



The second integral is constrained by the conditions 
K = p — p' , Ep > > + 4m2, and the constant 
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FIG. 1: Effective distortion distance D as a function of the 
isotropic oscillation-free coefficient — c'^"' in units of GeV~® 



for three neutrino energies: 17 GeV (solid line), 
(dashed line), and 100 TeV (dotted line). 
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TABLE XV: Estimated values of isotropic oscillation-frcc co- 
efficients at which the distortion distance D{E) for Cerenkov 
pair emission crosses 1000 km for energies 17 GeV, 43 GeV, 
and 100 TeV. Coefficient units are GeV'*"''. 



C is given by C = 2G|,M|(1 - 4s^ + 8s^)/{2Trf. In this 
second expression we have used the (5- function to perform 
one integration, introducing the convenient combinations 
K** = fc** + k''^ and k''^ = k'^ — fc'** with spatial spherical 
polar angles {0,(f>) and {9',4)'). 

The integral (157) provides the energy loss per distance 
in the presence of Dirac-type Lorentz- violating operators 
of arbitrary mass dimension, including anisotropic ef- 
fects. To gain some feeling for this integral and to extract 
estimated bounds on Lorentz violation from Cerenkov 
spectral distortion, consider the isotropic oscillation-free 
model with neutrino energy given by Eq. (86) in the 
preferred frame. For this special case, Eq. (155) shows 
the factor q/q^ becomes (l,p), so within our assumptions 
all the Lorentz violation lies in the kinematics. Given 
values of the coefficients a^'^\ c^'^^ and a specific neu- 
trino energy E^, the integral (157) can be numerically 
evaluated and used to calculate the effective distortion 
distance, defined by 



D{E) 



E 



dE/dx 



(158) 



For fixed energy, this distance is a rapidly falling function 
of the coefficients cS'^^ and c^*^^ , as illustrated in Fig. 1 for 
the case d = 10. It therefore offers a useful measure of 
the spectral distortion caused by Lorentz violation. 

No substantial spectral distortion is observed at the 
mean OPERA energy of 17 GeV or ranging up to ener- 
gies of over 40 GeV [14] , and no depletion of atmospheric 



neutrinos or antincutrinos is detected in IceCube for en- 
ergies over 100 TeV [132]. Requiring that the distortion 
distance D at these energies is 1000 km or greater, we 
can extract conservative limits on isotropic oscillation- 
free coefficients for Lorentz violation. Table XV shows 
estimates for the bounds obtained by taking one coeffi- 
cient to be nonzero at a time. The bounds are one sided 
except for the 100 TeV results for a^^^, a^^), and a^^\ 
which are on the modulus of the coefficients. The val- 
ues include ones substantially tighter than the results for 
the coefficients found in time-of-fiight experiments, which 
arc given in Table XIII. Note that a complete calculation 
at 100 TeV is likely to produce sharper limits because at 
these energies pair emission is expected to be only a small 
contribution. For example, only about 3% of the on-shell 
emission generates electron-positron pairs. 
An explicit bound on a coefficient for Lorentz viola- 
tion has recently been obtained from an analysis of the 
decay rate and energy loss for Cerenkov pair emission 
[79]. This work assumes that Lorentz violation arises 
only from the isotropic oscillation-free coefficient c'-'*^ in 
the minimal SME, leading to an isotropic constant shift 
5v = —cS^^ of the muon-ncutrino group and phase veloc- 
ities. The analysis finds the best limit on this type of 
Lorentz violation comes from IceCube measurements of 
100 TeV neutrinos, translating to the one-sided bound 
c'-^' > —8.5 X 10^^^. Also within this particular one- 
coefficient model, the ICARUS collaboration reports a 
limit based on the nonobservation of Cerenkov-emission 
products in the same neutrino beam [135], which corre- 
sponds to the one-sided bomid > —2 x 10^**. Both 
these limits are many orders of magnitude tighter than 
the nonzero negative value of c^^^ implied by the OPERA 
result for this one-coefficient model and listed in the first 
row of Table XIII. 



VIII. SUMMARY 

In this work, we study the effects of Lorentz and CPT 
violation on the behavior of fermions, with emphasis on 
neutrinos. The starting point in Sec. II is the construc- 
tion of the general quadratic Lagrange density (3) for the 
propagation and mixing of N species of fermions. This 
permits Lorentz- violating operators of arbitary mass di- 
mension to be classified and enumerated. A procedure to 
obtain the leading-order hamiltonian (19) for fermions is 
described. 

In Sec. Ill, the general theory is specialized to extract 
the Weyl hamiltonian (32) describing the propagation 
and mixing of three flavors of left-handed neutrinos in the 
presence of Lorentz- and CPT- violation involving opera- 
tors of arbitrary mass dimension. Block diagonalization 
of the Weyl hamiltonian at leading order in mass and 
Lorentz violation generates the effective hamiltonian (51) 
describing dominant modifications to neutrino propaga- 
tion and mixing. This key result depends on 10 sets of 
coefficients for Lorentz violation that enter the Lorentz- 
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TABLE XVI: Summary of spherical coefficients. 



32 

violating piece (44) of the effective hamiltonian. These 10 
sets and some of their features are listed in Table I. Six 
involve Dirac-type operators for propagation and mix- 
ing, while the other four are Majorana type and describe 
neutrino-antineutrino mixing. 

Since violations of rotation symmetry are a central fea- 
ture of many searches for Lorentz violation, it is useful to 
perform a decomposition of the effective hamiltonian and 
the coefficients using spherical harmonics. The resulting 
10 sets of spherical coefficients for Lorentz violation are 
presented in Sec. IV. Six are of Dirac type and four of 
Majorana type. This analysis also reveals that the fun- 
damental experimental observables for Lorentz violation 
in the neutrino sector comprise only four sets of effective 
spherical coefficients, built from the 10 basic sets accord- 
ing to Eq. (57). The properties of all 14 coefficient sets 
are given in Table II. 

Various special theoretical scenarios for the spherical 
coefficients can be countenanced, leading to different ex- 
perimental predictions. Section V presents several limit- 
ing cases of the general formalism. We begin in Sec. V A 
by matching to the renormalizable sector of the SME, 
revealing some qualitatively new effects that appear at 
leading order in both mass and Lorentz violation. Some 
properties of the renormalizable coefficients are given in 
Table III. 

In Sec. VB the limit of massless models is described. 
This class of models is of interest in part because it of- 
fers the potential for an alternative description of neu- 
trino oscillations without invoking mass. The massless 
coefficients and some properties are given in Table IV. 

Another scenario of potential interest, discussed in Sec. 
V C, is flavor-blind Lorentz violation. In these models the 
different neutrino flavors are assumed to experience the 
same effects, which is a reasonable approximation under 
some experimental conditions. For example, the limit 
of oscillation-free propagation is of relevance for certain 
types of searches for Lorentz violation, such as time-of- 
flight experiments. Table V lists the flavor-blind and 
oscillation-free coefficients and some of their features. 

In Sec. VD we consider diagonalizable models, in 
which all terms in the effective hamiltonian are assumed 
to be simultaneously diagonalizable. These models of- 
fer comparatively simple options for model building with 
Lorentz violation, while avoiding the complications that 
appear for general neutrino mixings. The diagonalizable 
coefficients are tabulated in Table VI. 

Finally, various types of isotropic models are studied 
in Sec. VE. Isotropy can be enforced only in a pre- 
ferred inertial frame, with anisotropies appearing in other 
boosted frames. These models are particularly simple be- 
cause they cannot have Majorana mixings and their Dirac 
terms must be isotropic. We discuss generic isotropic 
models and their restriction to isotropic diagonalizable 
models and to isotropic oscillation-free models. The co- 
efficients for each case are shown in Table VII. 

Table XVI summarizes the various spherical coeffi- 
cients for Lorentz violation studied in this work. The 
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first column of this table lists the theoretieal scenario, 
the second lists the relevant sets of spherical coefficients, 
the third provides the range of operator mass dimension 
d allowed for each coefficient set, and the last column 
indicates the number of independent coefficients for each 
value of d. The first four rows concern the effective spher- 
ical coefficients that are the fundamental observables for 
any experiment in this framework. The 10 following rows 
list the 10 basic sets of spherical coefficients arising in the 
general formalism, which separate into six sets of Dirac 
type and four sets of Majorana type. The rest of the 
table concerns the various limiting theoretical scenarios 
discussed in the text. 

In the penultimate sections of the paper, we study ex- 
perimental implications of the theoretical framework and 
use existing data on neutrino oscillations and propaga- 
tion to extract constraints on coefficients for Lorentz and 
CPT violation. Section VI treats effects on mixing. Two 
experimentally relevant limits are analyzed. The first is 
the short-baseline approximation, for which oscillation ef- 
fects from all sources are small. General expressions for 
short-baseline oscillation probabilities are given in Eq. 
(98). These expressions are applied to results from the 
short-baseline experiments LSND and MiniBooNE to ex- 
tract maximal attained sensitivities to effects from flavor- 
mixing Lorentz-violating operators for d < 10. The re- 
sults for generic spherical coefficients for Lorentz viola- 
tion are listed in Table VIII, while those for isotropic co- 
efficients are given in Table IX. The second limit is the 
perturbative approximation, in which the baseline can be 
arbitrary but the coefficients for Lorentz and CPT vio- 
lation are assumed small compared to the 3i^SM masses. 
We derive the mixing probabilities up to second order 
in Lorentz violation, and we present methods to ana- 
lyze data from long-baseline experiments. Numerical val- 
ues of beam-dependent factors relevant to the MINOS, 



OPERA, and T2K experiments are tabulated in Table X. 
As a simple example, we consider the limit of two-flavor 
mixing and discuss some asymmetries offering sensitivity 
to CPT violation. 

In Sec. VII, we discuss several types of kinematic tests 
that are independent of mixing. Effects in these tests are 
controlled by oscillation-free coefficients for Lorentz vio- 
lation. We consider various time-of-flight measurements, 
including the OPERA and MINOS experiments, earlier 
studies at Fermilab, and the supernova SN1987A. These 
are used to extract constraints on generic oscillation- 
free coefficients for d < 6 and on isotropic oscillation- 
free coefficients for d < 10. The results are collected 
in Tables XI, XII, and XIII. We also obtain threshold 
constraints from the observation of high-energy neutri- 
nos by IceCube, presented in Table XIV, and we derive 
the estimated order-of-magnitude bounds from neutrino 
Cerenkov radiation shown in Table XV. 

The analysis in this paper provides a general theoret- 
ical framework for the treatment of neutrino propaga- 
tion and mixing, along with a guide to its application in 
laboratory experiments and astrophysical observations. 
We see that searches involving neutrino propagation and 
mixing offer excellent sensitivities to numerous distinct 
types of Lorentz and CPT violation. Many coefficients 
for Lorentz violation remain unconstrained, so a substan- 
tial region of untouched territory is open for future ex- 
ploration using laboratory and astrophysical techniques. 
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